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SUMMARY 


Development of a computational method for prediction of external store carriage 
characteristics at transonic speeds is described. The geometric flexibility required 
for treatment of isolated and underwing, pylon-mounted stores is achieved by comput- 
ing finite difference solutions on a five-level embedded grid arrangement. A 
completely automated grid generation procedure facilitates applications. Store 
modelling capability consists of bodies of revolution with multiple fore and aft 
fins. A body-conforming grid improves the accuracy of the computed store body flow 
field. A nonlinear relaxation scheme developed specifically for modified transonic 
small disturbance flow equations enhances the method's numerical stability and 
accuracy. As a result, treatment of lower aspect ratio, more highly swept and 
tapered wing planforms is possible. A limited supersonic freestream capability is 
also provided. Pressure,' load distribution, and force/moment correlations show good 
agreement with experimental data for several test cases. A detailed computer program 
description for the Transonic Store Carriage Loads Prediction (TSCLP) Code is 
included in Appendix D. 
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INTRODUCTION 


Prediction of external store carriage characteristics at transonic speeds re- 
quires computations for rather complex geometries. Wing, fuselage, pylon, and store 
body and fin components each need to be modelled. While methods to obtain full 
potential, Euler, and Navier Stokes solutions for relatively simple geometries are 
maturing at a rapid pace, transonic small disturbance (TSD) formulations are still a 
practical alternative for treatment of these more complex configurations. 

The NASA/Grumman Transonic Wing-Body Code (Refs. 1,2,3) represents the 
state-of-the-art for reliable TSD analysis of complex aircraft. An attempt to extend 
similar wing/fuselage methodology to treat wing/fuselage/pylon/store geometries 
(Refs. 4,5) attributed poor isolated body normal force correlations to the TSD 
formulation. Since approaches emphasizing the use of more exact flow equations 
(Refs. 6,7,8) are difficult to implement and require further development for 
practical three-dimensional applications, a more accurate TSD formulation was 
developed for treatment of store body shapes (Ref. 9). This was accomplished by 
solving TSD flow equations on grids which conform to the store body shape, subject to 
exact (inviscid) store body surface boundary conditions. 

For this effort, refined TSD approaches (Refs. 1,2, 3, 9) have been combined for 
treatment of isolated and underwing, pylon-mounted stores. In the resulting Trans- 
sonic Store Carriage Loads Prediction (TSCLP) code, geometric flexibility is achieved 
by computing solutions on a five-level embedded grid arrangement. In addition, a 
nonlinear finite difference relaxation scheme developed specifically for modified TSD 
flow equations enhances numerical stability and accuracy. As a result, treatment of 
lower aspect ratio, more highly swept and tapered wing planforms is possible. A lim- 
ited supersonic freestream capability is also provided. 


3 
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NOMENCLATURE 


a,u,v 

A 

b/2 

c 


ave 


'f, ave 


m 

C 

m 

CMX 

CMY 

CMZ 

C n 

C N 


f 

f 

F 

i » J » k 


Coefficients in upwind difference formulae 
Area 

Semi span 
Chord 

Wing average chord (reference area 1 span) 
Sectional drag coefficient 
Drag coefficient 

Sectional skin friction coefficient 

Average turbulent skin friction coefficient 

Skin friction coefficient 

Sectional lift coefficient 

Rolling moment coefficient 

Lift coefficient 

Sectional moment coefficient 

Pitching moment coefficient 

Moment coefficient about x-axis 

Moment coefficient about y-axis 

Moment coefficient about z-axis 

Yawing moment coefficient 

Normal force coefficient 

2 

Pressure coefficient, (p-p )/ 

Store body cross-section axial force coefficient 
Axial force coefficient 

Store body cross-section side force coefficient 
Side force coefficient 

Store body cross-section normal force coefficient 

Normal force coefficient 

Airfoil section shape function 

Store body cross-section force coefficient 

Force Coefficient 

Grid indices 
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L 

M 

M 

n 


n x’ n r’ n 8 

W n z 

R 

Re 

RMAX 


s 

S 


T.U.V 

u ,v ,w 

oo oo 


oo 


V 

V 

oo 

X 


x,r,0 


x,y,z 

z 

a 


& 


s 


Y 

r 

6 


fin 


A 

C 

T| 

A 


i,f\ ,0 


Unit vectors 
Reference length 

Moment coefficient about reference center 
Freestream Mach number 
Surface normal 

Cylindrical components of surface normal 

Cartesian components of surface normal 

Store body radius 

Reynolds number 

Store body maximum radius 

Arc length 

Surface area 

Coefficients of governing flow equation 

Cylindrical components of freestream velocity 

Local velocity 

Freestream velocity 

Vector coordinate 

Cylindrical coordinates 

Cartesian coordinates 

x,r location in physical domain of conformal mapping 

Aircraft or isolated store angl e-of-attack relative to freestream 

Store pitch angle relative to aircraft (positive, nose up) 

Pylon yaw angle relative to aircraft (positive, leading edge 
outboard) 

Store yaw angle relative to aircraft (positive, nose outboard) 

Specific heat ratio 
Ci rcul ati on 

Store fin deflection angle (positive, leading edge counter-clockwise 
looking upstream) 

Incremental quantity, or mesh cell size 

£,ti location in transformed domain of conformal mapping 

Wing or fin non-dimensional spanwise location (2y/b or 2r/b) 

Local sweep angle of constant percent chord line 
Vi scosity 

Computational coordinates in Cartesian grids 
Computational coordi nates i n cylindrical grids 
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p 


0 

W 


>0 

nn ,y ss 
s 


oo 


Density 

Perturbation potential, or store pressure-tap roll orientation angl 
Central and upwind second derivatives of perturbation potential 
Store roll rate, p/V (radians per unit length) 

oo 

Infinity 


Subscri pts 


b,body 

c 

f ,fin 

LE,TE 

p, pylon 

REF 

REFS 

s, store 

u,l 

w, wi ng 
wet 

x, r,0 
x,y,z 


Store body 

Viscous crossflow quantity 
Store fin 

Leading edge or trailing edge 
Pylon 

Configuration or isolated store reference parameter 

Store reference parameter 

Store 

Upper or lower 
Wing 

Wetted area 

Partial derivatives in cylindrical coordinates 
Partial derivatives in Cartesian coordinates 


Superscripts 

D Dummy val ue 

+ New or updated value 
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COMPUTATIONAL METHOD 


Inputs to the TSCLP code consist of configuration geometry, freestream flow 
conditions, and number of solution iteration cycles. No additional solution or grid 
generation parameters are required. The grid generation procedure has been 
completely automated to facilitate applications. Treatment of wing/fuselage/pylon 
geometry is similar to that found in the basic NASA/Grumman Transonic Wing-Body Code 
(Refs. 1, 2), with enhancements as noted. Modelling capability for isolated and 
underwing, pylon-mounted stores is described below. 

An example of the complexity of store geometry that can be modelled is the GBU- 
15-CWW, shown in Fig. 1. In general, the store body must be axisymmetric and may be 
either sharp or blunt nosed. Two sets of fins may be input, each set consisting of 

i 

from one to four identical fins at different angular locations. When two sets of 
fins are input, they must be ’Of the fore and aft type, and they must have the same 
number of fins at the same angular locations. Fore and aft fins can have different 
planforms and (symmetric) airfoil sections. Individual, all-moveable fin deflection 
and quasi-steady, isolated store roll rate capabilities have also been incorporated. 



Transonic Small Disturbance Formulation 

Wing/fuselage/pylon (Cartesian grid) calculations employ the following 
"modified" TSD flow equation (Refs. 1,2): 

Td> + Ud> + Vd> + <J> = 0 (1) 

y xx y xy v yy y zz 1 ' 


9 




where 


T = 1 - M 2 - (y + 1) M 2 4> - M 2 <J> 2 

00 v ' 7 00 Y x Z 00 Y x 

U = -2M 2 <b, 

® r y 

V - 1 - (y - 1) M 2 * x 

Store body and fin (cylindrical grid) calculations employ a similar equation (Refs. 
3, 9): 

T d> + Ud> + V(b + ~ d> + — r d> Q n = 0 ^ 2 ^ 

y xx y xr y rr r y r r 1 v 00 

where 

T = 1 - M 2 - (y + 1)M 2 <J> - ^4-^ M 2 <j> 2 

00 v ' 9 00 Y x 2 00 Y x 

U = -2M 2 (a sin0 + A ) 

oo v r r 

V - 1 - (Y - 1) H 2 * x 

These TSD flow equations are "modified" by the retention of terms neglected by 

2 

"classical" TSD flow equations (Ref. 10). The retention of the <f> x <|> term provides 

a better approximation to the transition between subsonic and supersonic flow. 

The <J> <b <j> d> and <(><(> <J> (J> terms are retained to more accurately resolve 

y xy , x yy r xr , x rr 

shock waves with appreciable sweep in the x-y (wing) and x-r (fin) planes, 
respectively. 

Lifting surfaces are treated using a small disturbance, planar approximation. 
For example, wing boundary conditions are imposed on the wing reference plane 
(z=z w .j n g) and, for an airfoil upper and/or lower surface section shape f(x), are 
given by: 

± <t> z = - " x (3) 

with the assumptions that: 
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n x = t [tan 1 (f x ) - a] 



In the wing fine grid, where wing leading and trailing edges are always located at 
grid mesh cell midpoints, it is possible to use a more accurate approximation for n x : 

n x ■ 1 < f x • “> 

Pressure coefficients on the wing surface are defined as: 

Cp - -2 -d - *1 (4) 

Pylon surfaces receive similar treatment. Boundary conditions are imposed on 
the pylon reference plane (y=ypyL0N^» ar| d are given by: 

±<j>y = - (cosa n x + sina n z ) (5) 

where 

n x - T [tan _1 (fx) - 6 p ] 

"y " * 1 

n = ± tan (A) tan *(f ) 

Z A 

Here the pylon surface normals are assumed to be perpendicular to swept, constant 
percent chord lines. Pressure coefficients on the pylon surface are defined as: 

Cp - -2* x - - (1 - »l) il (6) 

Body-type components require a different approach. The fuselage is treated 
using an approximate, constant cross-section boundary condition support surface 
(Ref. 1). For the store body, a more exact boundary condition is used: 
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K + *x) n x + ( V * + * r ) n r “ 0 


(7) 


where 

"x ‘ ' 37 / '/ 1 * © 2 

"r’ 1 /A* 

n 9 ■ 0 

and 

= cos (a + a $ ) cos (B $ ) 

V m = sin ($ s ) cos (0) + sin (a + a $ ) cos (0 g ) sin (0) 

- sin (a + a $ ) cos (0 $ ) cos (0) - sin (0 ) sin (0) 

This boundary condition is consistent with the body-conforming grids used for 
modelling the store body shape. The isentropic form of the pressure coefficient is 
also used on the store body surface: 

2 f r y - 1 2 2 n T 7 (T -1 ) 

C P V 1 ♦ 3 L F iH . <* - v >i -y (») 

where 

v 2 - (u„ * ♦/ * (V. * V ) 2 * w. * ? + s > 2 

Store fin surfaces are treated using a small disturbance, planar approximation 
similar to that for wings and pylons. Fin boundary conditions are imposed on the 
fin reference plane (0 = ©pj^)* and are given by: 

± £ <fV> = “ ( u n v + v nj t (W - r w ) (9) 

T o v oo X © T oo S 7 v ' 


where 
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n x - t [tan-‘(f x ) - «„„] 
n r • tan(A) tan -1 (f x ) 
n 9 - ± 1 

and pressure coefficients on store fin surfaces are defined as: 

c p - -2+ x - t r 2 - (1 - Mf) <>l (10) 

Grid System Arrangement 

A five-level, Cartesian/cylindrical, embedded grid arrangement is employed to 
facilitate treatment of wing/fuselage/pylon/store configurations. Geometry input 
verification plots and grid system arrangements are shown for the isolated GBU-15- 
CWW store in Figs. 2 and 3, for a Douglas wing/pylon/store test configuration in 
Figs. 4 and 5, and for a Nielsen wing/fuselage/pylon/store test configuration in 
Figs. 6 and 7. These figures supplement the general discussion of grid system 
arrangement which follows. Adetailed description of the automated grid generation 
procedure can be found in Appendix A. 

No provision is made for representing wing or fuselage surfaces in cylindrical 
grid systems, or for representing store body and fin surfaces in Cartesian grid 
systems. The pylon surface is modelled in both types of grids. Current grid setup 
will treat wing/store gaps as small as one store diameter, excluding fins, or one- 
half store radius, including fins, whichever is larger. Smaller gaps may compromise 
the grid generation procedure or result in erroneous solutions. 

First, a Cartesian global coarse grid is arranged about the entire config- 
uration. This grid is used to impose far field boundary conditions and to compute a 
coarse wing/fuselage/pylon ( i . e . , store off) flow field, which in turn provides a 
starting solution for embedded grid systems. A global coarse grid inner boundary 
surrounding the store is then defined. Inside this boundary no further global 
coarse grid computations are performed. 
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Figure 4 Input Geometry Verification Plots for Douglas Wing/Pylon/Store Configuration 
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Figure 5 Grid System Arrangement for Douglas Wing/Pylon/Store Configuration 
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Figure 6 Input Geometry Verification Plots for Nielsen Wing/Fuselage/Pylon/Store Configuration 
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Figure 7 Grid System Arrangement for Nielsen Wing/Fuselage/Pylon/Store Configuration 











Next, a Cartesian Wing/Store Interaction (WSI) grid is embedded about the 
store, within the global coarse grid. This is a medium density grid which functions 
primarily as a means of communication between Cartesian and cylindrical grid 
systems. An overlap region is created between the global coarse grid inner boundary 
and the WSI grid outer boundary, so that each can be updated using flow field 
information from the other grid. A WSI grid inner boundary surrounds the store 
(more closely than the global coarse grid inner boundary). Inside this boundary no 
WSI grid computations are performed. 

A coarse cylindrical C-grid is then embedded about the store body, within the 
WSI grid. The use of a body-conforming grid improves the accuracy of the computed 
store body flow field. A conformal mapping provides the basic C-grid transformation 
(Ref. 11): 

z « in (1 - cosh c). (11) 

Stretching and shearing transformations are also employed. An overlap region is 

f 

created between the WSI grid inner boundary and the coarse C-grid outer boundary, so 
that each may be updated using information from the other grid. The C-grid is 
further modified so that constant coordinate grid lines approximate the fin tip 
vortex streamline locations. This facilitates treatment of fins and their wakes. 

As in the basic Transonic Wing-Body Code, an embedded wing fine grid system is 
used to improve solution accuracy near the wing. In a similar fashion, a fine 
cylindrical C-grid is embedded within the coarse C-grid, to improve the accuracy of 
the computed flow field near the store. 

For subsonic freestreams (M^ < 1.0), transformations locate outermost grid 
boundaries at infinity, where appropriate far field boundary conditions are 
applied. For supersonic freestreams (M m 1.0), outermost grid boundaries are 
located a finite distance from the configuration, where supersonic inflow, outflow, 
and radiation-type boundary conditions are employed. 

Solution Algorithm 

The following discussion gives a general overview of the solution algorithm. 

A detailed description of the finite difference approximations employed can be found 
in Appendix B. 
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At each grid point, finite difference approximations are substituted for terms 
appearing in the governing flow equations, Eqs. 1 and 2. An upwind, rotated 
difference scheme (Ref. 12) provides the proper domain of dependence at supersonic 
points. A variation of this scheme, developed specifically for modified TSD flow 
equations (Ref. 3), determines the rotation from the coefficients T,U,V rather than 
from local flow angularities. This greatly enhances the method's numerical 
stability and accuracy. 

As a result, treatment of lower aspect ratio, more highly swept and tapered 
wings is possible. In conjunction with appropriate inflow, outflow, and radiation- 
type boundary conditions (Ref. 13), a supersonic freestream capability is also 
provided. This works well for wing/fuselage combinations, but only limited success 
was achieved for supersonic treatment of stores. 

Communication across embedded grid boundaries is accomplished via Neumann-type 
boundary conditions for central first derivatives and, where required, upwind second 
derivatives of the potential in a direction normal to the computational grid 
boundary. An overlap region is always created between embedded grid boundaries, so 
that each can be updated using flow field information from the other grid. The 
desired flow field quantity is obtained from the other grid by linear interpolation 
between mesh cell corner points. Neumann-type boundary conditions were found to 
improve code covergence and minimize shock reflections at embedded grid boundaries, 
relative to Dirichlet-type boundary conditions. 

Treatment of wing, fuselage, and pylon surfaces is similar to that found in the 
basic Transonic Wing-Body Code (Refs. 1,2), with the following exceptions. To 
retain second order accuracy, pylon surface boundary conditions are imposed in the 
wing fine grid using a Z-scheme (Ref. 14). To maintain numerical stability at the 
wing/pylon junction, second order accurate wing and pylon junction boundary 
conditions are replaced in all grids by a combination of first and second order 
accurate formulae. 

The accuracy of the computed store body flow field is improved by the use of 
body-conforming grids. First order accurate boundary conditions are used to set 
potentials on the store body surface. Fin surfaces are treated using a small 
disturbance, planar approximation similar to that used for the wing and pylon. At a 
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store body/fin or store body/pylon junction, the second order accurate fin or pylon 
surface boundary condition is replaced by a first order accurate formula. This 
formula is combined with the store body surface boundary condition, so that a single 
value of the potential can be set in the junction. 

The embedded grid solution process is divided into three phases. Computations 
during each phase proceed as follows. 

A "coarse grid" solution is first obtained for wing/fuselage/pylon ( i . e . , store 
off) geometry in the global coarse grid. This initial phase allows for rapid 
propagation of disturbances and also provides starting.flow fields for the other 
grid systems. Each iteration cycle consists of a global coarse grid relaxation 
sweep and, for supersonic freestreams, an outer boundary update. For isolated 
stores, global coarse grid potentials are merely set to zero. 

Next, an "intermediate grid" solution is obtained for the complete wing/ 
fuselage/pylon/store geometry. This phase allows for rapid calculation of the store 

f 

body and fin flow field as well as any airframe/store interference effects. It also 
provides starting flow fields for the wing and store fine grid systems. First, WSI 
grid and coarse C-grid potentials are initialized based on the "coarse grid" 
solution. Each "intermediate grid" iteration cycle then consists of inner and outer 
grid boundary updates and grid relaxation sweeps for the global coarse grid, WSI 
grid, and coarse C-grid. After each WSI grid relaxation sweep any global coarse 
grid wing surface points coincident with the global coarse grid inner boundary are 
set based on WSI grid values. Potentials at these points are then held fixed during 
each global coarse grid relaxation sweep. 

The third phase consists of the "fine grid" solution process. First, the wing 
fine grid and fine C-grid are initialized based on the "intermediate grid" 
solution. Each "fine grid" iteration cycle then includes inner and outer grid 
boundary updates and grid relaxation sweeps for all grid systems. After each wing 
fine grid relaxation sweep all global coarse grid and WSI grid wing surface points 
are set based on wing fine grid values. Potentials at these points are then held 
fixed during each global coarse grid and WSI grid relaxation sweep. After each fine 
C-grid relaxation sweep all coarse C-grid store body, store fin, and pylon surface 
points are set based on fine C-grid values. Potentials at these points are then 
held fixed during each coarse C-grid relaxation sweep. 
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Force and Moment Coefficients, Including Viscous Effects 

After a solution has been obtained, surface pressures are integrated to yield 
load distributions and force and moment coefficients. For cases where the wing 
finite difference boundary layer calculation is not activated, estimated wing upper 
and lower surface section skin friction coefficients are obtained in the same manner 
as for the fuselage, using the Prandtl-Schlichting formula (Ref. 15) corrected for 
compressibility effects: 

c f,ave = (1 + °* 028 M J °- 455/ [ 1 og (Re)] 2 ' 58 (12) 

where the Reynolds number is based on local wing chord. 

In general, the pylon surface will not be represented in its entirety in any 
one grid system. It therefore becomes necessary to piece together pylon coefficient 
contributions from several grid systems. First, contributions from that portion of 
the pylon located in the fine C-grid are computed. Next, contributions from that 
portion of the pylon surface located in the wing fine grid (but not in the fine 
C-grid) are considered. Similarly, contributions from portions of the pylon surface 
in the global coarse grid and MSI grid are considered, as required, until the entire 
pylon surface has been accounted for. Pylon section skin friction coefficients are 
also estimated using Eq. 12 and a Reynolds number based on local pylon chord. 

Store body loads are calculated based on computed inviscid pressure 
coefficients. A skin friction coefficient for the store body, based on wetted area, 
is estimated using Eq. 12 and a Reynolds number based on store body length. Fin 
section skin friction coefficients are estimated using Eq. 12 and a Reynolds number 
based on local fin chord. Store body viscous crossflow effects (Refs. 16, 17) are 
also estimated, based on crossflow Reynolds number, flow angularity with respect to 
the freestream, and store body fineness ratio. This viscous crossflow estimate is 
strictly valid for isolated store bodies only, since it does not account for fin 
and/or airframe interference effects. 

A detailed description of the computation of force and moment coefficients can 
be found in Appendix C. 
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EVALUATION OF RESULTS 


Results computed by the TSCLP code for a variety of test cases are evaluated by 
comparisons with more exact methods and with experimental data. Test cases con- 
sidered include isolated stores and simple wing/fuselage combinations, as well as 
complete wing/fuselage/pylon/store configurations. Component interference effects 
are also examined. 

Isolated Stores 

Test cases for isolated stores include several body-alone geometries and one 
which features multiple fore and aft fins. All calculations were made using 200 "me- 
dium grid" iterations and 200 "fine grid" iterations. 

The first two test cases are the National Transonic Facility (NTF) 5° Cali- 

* 

bration Cone at M =0.6 and a=0°, and the Pathfinder I Nosecone at M =0.84 and a=0°. 

OO • * oo 

Body shapes appear in Fig. 8. Note the rather sharp nose of the former, and the rel- 
atively blunt nose of the latter. The figure compares calculated body pressure dis- 
tributions with results (Ref. 9) obtained using a two-dimensional (axi symmetric) full 
potential method. Results compare well for both test cases. 

The next test case is the isolated Nielsen generic store at M ro =0.925 and a=5°. ■ 
This store is shown in the carriage position in Fig. 6. Calculated body pressure and 
normal force load distributions for the isolated store are compared to experimental 
data (Ref. 18) in Fig. 9. Good correlation with data is observed. Note the small, 
almost imperceptible change in the load distribution due to estimated viscous cross- 
flow effects. 

Total integrated normal force and pitching moment coefficients are presented in 
Table I. Computed inviscid coefficients agree well with slender body theory. Upon 
integration, estimated viscous crossflow effects are significant and markedly improve 
the correlation with experiment. 
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Figure 8 Comparison of Isolated Body Pressure Predictions with Axisymmetric Full Potential 
Calculations 



Figure 9 Isolated Nielsen Store: Pressure and Load Distribution Correlations, M a = 0.925, a = 5° 
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Table I Isolated Nielsen Store: Force and Moment Correlation at Mg, = 0.925 and a = 5° 



Slender 

Body 

Theory 

TSCLP 

(Inviscid) 

TSCLP 

w/Viscous 

Crossflow 

Experiment 

C N 

0.175 

0.176 

0.239 

0.221 


0.582 

0.578 

0.558 

0.462 


6096-026 

Attempts to analyze stores at supersonic speeds met with limited success. Lack 
of convergence for most cases is attributed to the C-grid relaxation scheme. In the 
forward portion of this grid, sweeping from the body surface to the grid outer 
boundary is commensurate to marching in an upstream direction. To compound the 
situation, grid stretching in this region is rather severe. Further development is 
required for reliable treatment of stores at supersonic speeds. 

f 

A converged solution was, however, obtained for a NACA research body at M =1.0 

t oo 

and a=8°. The body shape appears in Fig. 10, which compares calculations at M^=0.99 
and Mj=1.00 to experimental data (Ref. 19) at M m = 1.0. Corresponding body normal 
force correlations appear in Table II. Predictions agree fairly well with data. 
Despite the difficulties mentioned earlier, these calculations demonstrate that the 
supersonic inflow, outflow, and radiation-type boundary conditions implemented are 
i ndeed vi abl e . 


Table II NACA RM L53H04 100-inch Body: Normal Force Correlation at M^, s 1.0 and a = 8° 



Slender 

Body 

Theory 

TSCLP 

(Inviscid) 

Moo = 0.99/1.00 

TSCLP 
w/Viscous 
Crossflow 
Mqo = 0.99/1.00 

Experiment 

c N 

0.074 

0.073/0.075 

0.160/0.162 

0.190 


6096-027 


The final isolated store test case is the GBU-15-CWW store featuring multiple 
fore and aft fins. Calculations were made at M =0.95 and a=6° for several fin ar- 

oo 

rangements. Calculated pressure distributions are compared to experimental data 
(Ref. 20) in the following figures: 
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Fig. 11: Body pressures, wings and canards off 
Fig. 12: Body pressures, wings and canards on 
Fig. 13: Wing pressures, canards off 
Fig. 14: Wing pressures, canards on 
Fig. 15: Canard pressures. 

Overall, correlation with data is very good. A forebody double shock system which is 
predicted for the fins off arrangement is not present in the data. This discrepancy 
affects the canard pressure comparisons as well, although the latter are fairly good 
considering the sparse number of grid points used to represent the canard surfaces 
(see Fig. 3). The effect of the fins on the body is accurately predicted,, as is the 
downwash effect of the canards on the wings. 
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Figure 11 GBU-15-CWW Store: Body Pressure Distribution Correlations, Wings and Canards Off, 
Moo = 0.95, a = 6° 



Figure 12 GBU-15-CWW Store: Body Pressure Distribution Correlations, Wings and Canards On, 
Mqo = 0.95, a = 6° 


Total integrated normal force and pitching moment coefficients for the 
GBU-15-CWW store are compared with data in Fig. 16. Results are shown as a buildup 
of body, wing, and canard components, and include body viscous crossflow effects. 

The underprediction of wing-on normal force is attributed to sparse wing grid point 
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Figure 13 GBU-15-CWW Store: Wing Pressure Distribution Correlations, Canards Off, = 0.95, 
a = 6 ° 


distribution and to a not yet fully converged solution (complete convergence requires 
several times the number of iteration cycles typically used for the present, 
engineering calculations), as well as to flow phenomena not modelled by the governing 
flow equations. Canard normal force and pitching moment increments are again 
indicative of their downwash effect on the wings, and are properly predicted. 
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Figure 14 GBU-15-CWW Store: Wing Pressure Distribution Correlations, Canards On, M m = 0.95, 
a = 6° 


Store Carriage Configurations 

Results for two store carriage configurations are evaluated by comparisons with 
experimental data. All calculations were made using 200 "coarse grid" iterations, 
200 "medium grid" iterations, and 200 "fine grid" iterations. Since no boundary 
layer computations are provided for pylon and store surfaces, those for the wing were 
not employed. Store characteristics do, however, include the simple estimate of body 
viscous crossflow effects. 
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Figure 15 GBU-15-CWW Store: Canard Pressure Distribution Correlations, M ffl = 0.95, a = 6 V 


The first store carriage test case is the Douglas wing/pylon/store configuration 
shown in Fig. 4. Although no store surface pressure or loads data is available, wing 
surface pressure data obtained at a span station in the vicinity of the pylon (Ref. 
21) permits an evaluation of computed results. Correlations for the isolated wing 
and for the complete wing/pylon/store combination at M =0.75 and a=4° appear in Fig. 17 
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Figure 16 GBU-15-CWW Store: Force and Moment Correlation, M w = 
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Figure 17 Douglas Wing/Pylon/Store Configuration: Wing Pressure Distribution Correlations and 
Spanload Comparison, = 0.75, a = 4° 







Even for the isolated wing, correlation of pressures on the upper surface is 
poor. Airfoil section analyses using the more sophisticated 2-D GRUMFOIL method 
(Ref. 22) gave similar predictions, suggesting that test anomalies, rather than flow 
simulation inaccuracies, are to blame. Nevertheless, the predicted and measured 
effect of the pylon and store on wing upper surface pressures is a more pronounced 
leading edge expansion. 

Correlation with data in Fig. 17 is better on the wing lower surface. The pylon 
leading edge compression and subsequent expansion propagate onto the wing lower sur- 
face, as do the pylon trailing edge expansion and recompression. These effects occur 
in the presence of an overall acceleration through the converging/diverging region 
created by the wing and store. The very large expansion and strong shock predicted 
at the pylon trailing edge closure occur in the data to a lesser degree, possibly due 
to viscous effects and/or flow separation which have not been modelled. 

Computed wing spanload distributions are also shown in Fig. 17. The combined 
effect of pylon and store is a marked decrease in wing loading which is greatest at 
the wing/pylon junction itself. The side force carried on the pylon appears as a 
discontinuity in the wing spanload. 

The second store carriage test case considered is the Nielsen wing/fuselage/ 
pylon/store configuration shown in Fig. 6. Results were computed at M =0.925 and 

oo 

a=5° for isolated wing/fuselage and store components as well as for the wing/ 
fuselage/pylon and wing/fuselage/pylon/store combinations. Experimental data (Ref. 
18) is available for fuselage and store components only. 

Fuselage bottom centerline pressure distribution correlations appear in Fig. 18. 
For the wing/fuselage configuration, the wing lower surface pressure field propagates 
onto the fuselage bottom centerline as expected. With the pylon present this effect 
is enhanced slightly, and even more so with both pylon and store present. Pre- 
dictions correlate well with data. 

Store surface pressure correlations are presented in Fig. 19. Overall, agree- 
ment is good. While predicted expansions and compressions are slightly larger than 
those present in the data (as noted previously for the Douglas test configuration), 
the strong store/airframe interactions are properly predicted. 
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Figure 18 Nielsen Wing/Fuselage/Pylon/Store Configuration: Fuselage Bottom Centerline Pressure 
Distribution Correlations, M oo = 0.925, a = 5° 
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Figure 19 Nielsen Wing/Fuselage/Pylon/Store Configuration: Store Pressure Distribution 
Correlations, M a = 0.925, a = 5° 
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Pylon-mounted store axial load distribution correlations are shown in Fig. 20. 
Good correlation is shown for both normal force and side force load distributions, 
except for anomalies in predicted side force in the vicinity of the pylon and wing 
trailing edges. These anomalies are attributed to small misalignments of predicted 
inboard/outboard pylon and wing trailing edge shock locations (see Fig. 19). 
Unfortunately, relatively large loads result from rather small differences in surface 
pressures. Modelling of viscous effects and/or flow separation might correct this 
computational deficiency. 
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Figure 20 Nielsen Wing/Fuselage/Pylon/Store Configuration: Store Load Distribution Correlations, 
Mo, = 0.925, a = 5° 


Incremental forces and moments for the pylon-mounted store, relative to the 
isolated store, appear in Table III. While qualitatively correct, the predicted 
incremental effects of the airframe on the store lack the desired accuracy. Again, 
modelling of viscous effects and/or flow separation might improve predictions. 
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Table III Nielsen Wing/Fuselage/Pylon/Store Configuration: Correlation of Store Incremental 
Forces and Moments (Relative to Isolated Store) at Mq, = 0.925 and a s 5° 



TSCLP 

Experiment 

ac n 

-0.126 

-0.063 

AC m 

-0.510 

-0.712 

ACy 

0.158 

0.190 

AC n 

-0.694 

-0.294 


6096-028 


The Nielsen wing/fuselage/pylon/store configuration was tested at supersonic, as 
well as subsonic, Mach numbers. Although limited success was achieved for supersonic 
treatment of stores, calculations for this configuration did converge, possibly be- 
cause the underwing store is in a region of reduced, lower Mach number flow. Results 
computed at M^=l.l and a=5° are shown in Figs. 21-23. 
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Figure 21 Nielsen Wing/Fuselage/Pylon/Store Configuration: 
Fuselage Bottom Centerline Pressure Distribution 
Correlation, M m = 1.1, a = 5° 
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Figure 22 Nielsen Wing/Fuselage/Pylon/Store Configuration: Store Pressure Distribution 
Correlations, M tt = 1.1 9 a = 5° 



Figure 23 Nielsen Wing/Fuselage/Pylon/Store Configuration: Store Load Distribution Correlations, 
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The fuselage bottom centerline pressure distribution correlation is shown in 
Fig. 21. Predictions agree well with data. Store surface pressure comparisons ap- 
pear in Fig. 22. Strong shocks which are predicted near the back end of the store do 
not occur in the data. These discrepancies compromise the store axial load dis- 
tribution correlations, which appear in Fig. 23. Overall, correlation with data is 
fair, not nearly as good as for the subsonic (M^O.925) condition. 

Discrepancies in Figs. 21 through 23 are attributed to viscous effects and/or 
flow separation which have not been modelled, and also to the numerical deficiencies 
discussed earlier relative to supersonic treatment of isolated stores. In addition, 
it is possible that grid setup and embedded grid boundary treatment require further 
development to reliably predict store carriage characteristics at supersonic speeds. 

Low Aspect Ratio, Highly Swept and Tapered Wings at Supersonic Speeds 

Supersonic calculations were made for two low aspect ratio, highly swept and 

tapered wing geometries. Results illustrate the enhanced numerical stability of the 

upwind rotated finite difference scheme. Calculations were performed using 200 

# 

"coarse grid" iterations and 200 "fine grid" iterations. Wing boundary layer compu- 
tations were not employed. 

The first swept wing test case is the F-14 aircraft at M =1.3 and a=5°. The 

oo 

wing leading edge sweep is 68°. Computed pressure distributions are compared to 
experimental data (Ref. 23) in Fig. 24. Although some of the data appears erratic, 
several trends can be identified. Correlation with data is better at inboard 
stations than at outboard stations. Outboard, upper surface leading edge expansions 
are overpredicted. Also, a region of supercritical crossflow is observed in the data 
(terminated by the shock swept from the inboard leading edge location to the tip 
midchord location) but not at all in the calculation. Discrepancies are attributed 
to flow separation which is not modelled and to numerical deficiencies which are dis- 
cussed below. 

The second swept wing test case is the SC3 Demonstration Wing. Leading edge 
sweep is 65° inboard and 57° outboard. This wing was designed for supersonic 
maneuver using a full potential method. It features attached flow and controlled 
supercritical crossflow at the design conditions M =1.62 and a=12°. Computed 
spanwise pressure distributions are compared to experimental data (Ref. 24) in Fig. 

25. Again, the supercritical crossflow is not adequately resolved in the compu- 
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tations. This is attributed to transonic small disturbance flow equation limitations 
and/or to the need for higher spanwise grid density (there are only 18 grid span cuts 
from root to tip). The inability to resolve regions of supercritical crossflow 
should not pose problems when spanwise flow remains subsonic (e.g., for similar wing 
planforms at transonic speeds, or for less swept and tapered planforms at supersonic 
speeds). Even so, predictions shown in Fig. 25 are more accurate than panel method 
results (Ref. 25) which are also commonly employed at these flight conditions. 
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CONCLUDING REMARKS 


Correlations indicate that the TSCLP code provides reliable prediction of ex- 
ternal store carriage pressure/loading characteristics. Predictions for isolated 
stores are good, especially when coupled with simple estimates of store body loads 
due to viscous crossflow effects. For underwing, pylon-mounted stores, even simple 
2-D or axi symmetric modelling of viscous boundary layer and flow separation effects 
on wing, pylon, and store surfaces would improve loads prediction accuracy. Further 
studies are required to fully understand the limitations of the transonic small 
disturbance, embedded grid formulation. Beyond these limitations, the method will 
still provide a fundamental understanding of complex store/airframe interactions not 
available by other means. 

A finite difference relaxation scheme developed specifically for modified 
transonic small disturbance flow equations; allows treatment of lower aspect ratio, 
more highly swept and tapered wings, and provides a supersonic freestream capability 
(supersonic calculations for stores met with limited success). Numerical stability 
and accuracy is enhanced to such a degree that wing/fuselage solutions are readily 
obtained even for cases where the physics of the flow are well beyond the scope of 
the governing flow equation. As with any computational method, the user should 
assess the validity of a solution for each particular application. 
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APPENDIX A 

GRID GENERATION PROCEDURE 


Global Coarse Grid 

The Cartesian global coarse grid is arranged about the entire configuration. 
Initially, 71 grid points are evenly spaced in the x direction, over a finite region 
extending one half the total length upstream and downstream of the configuration: 


x - A q + Aj £ 


(13) 


with -1.0 £ £ $ +1.0. For subsonic freestreams, the first eight and the last eight 
mesh cells are redefined so as to stretch to upstream and downstream infinity: 


x 


= A q + Aj £ + A^ tan [ 



( 


(14) 


The constant A 3 is chosen so that next-to-last grid points are located twice the 
total length upstream and downstream of the configuration. 


For cases with a wing or fuselage, only the half space y i 0 is modelled: 

y “ tanh " 1 (n) (15) 

B 1 

with 0.0 £ p £ 1 . 0 , and 

z = Cj tan [\ 5 ) (16) 

with -1.0 £ £ £ +1.0. Total number of grid points is 26 in the y direction and 31 
in the z direction. The constant Bj is chosen so as to provide 18-1/2 mesh cells to 
the wing tip or, when no wing is present, 2 - 1/2 mesh cells to the fuselage side. 

The consant Cj is chosen to be the smaller of two values determined by the 
requirements that next-to-last grid points fall five wing average chord lengths 
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above and below the configuration, and, when a fuselage is present, that z mesh cell 
size at £=0 be equal to y mesh cell size at q»0. 


When a pylon is present, a small bump is added to the spanwise transformation 
from centerline to wing tip so as to ensure a global coarse grid line at the pylon 
spanwise location: 


y = -g- tanh -1 (n + An) . 


1 


(17) 


A combination of polynomial and exponential functions are used to define the bump: 
r / . , /e^^TIP - e^ri/n-rrp 2 

An = 64 [tanh(B r y pyL0N )- n pyL0N ]( R (18) 

e - 1 e - 1 

where 


YpYLON = ac l- ua ^ Pyl° n spanwise location 

n P Y|_oN = grid line desired at that location 

(closet grid line after applying Eq. 15) 

n T ip = tanh (B^ b W ING /2 ) 
and k is found by requiring that: 


e kn PYL0N /T1 TIP -1 1 



For isolated stores, the entire flow field is modelled: 

y - Bj tan ( J n) (19) 

with -1.0 £ n * +1.0. In the z direction Eq. 16 is used. Total number of grid 
points is 25 in both the y and z directions. The constants B^ and Cj are chosen so 
as to provide 2-1/2 mesh cells to the store body maximum radius. 

For supersonic freestreams, global coarse grid setup in the y and z directions 
is redefined, using cubic polynomials, so as to establish finite grid outer 
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boundaries in regions inboard of, outboard of, above, and below the configuration: 

y = b Q + bjH + b 2 n 2 + b 3 n 3 (20) 

z = c Q + CjC + c 2 C 2 + c^ 3 . (21) 

Global coarse grid setup near the configuration itself remains unchanged. 

Polynomial coefficients are determined by requiring continuous first and second 
order grid metrics where original and redefined regions abut, and by requiring a 
zero second order grid metric at the finite grid outer boundaries. 

When a store is present, the global coarse grid inner boundary is defined by 
those global coarse grid y and z grid lines which just enclose the store. This 
rectangular boundary is defined for the total upstream and downstream extent of the 
grid. 

Wing/Store Interaction Grid 

The Cartesian Wing/Store Interaction (WSI) grid is arranged about the store. 
Initially, 51 grid points are evenly spaced in the x direction, over a finite region 
extending one half the total store length upstream and downstream of the store nose 
and tail. Eq. 13 is used with -3.0 £ £ £ +3.0. For subsonic freestreams, the first 
eight and the last eight mesh cells are redefined using Eq. 14 so as to stretch to 
upstream and downstream infinity. The constants A3 are chosen so that next-to-last 
WSI grid points coincide with next-to-last global coarse grid points. For super- 
sonic freestreams, a finite transformation is employed: 

x = A q + A x ? + A 3 (? - ^) 3 . (22) 

The constants A3 are chosen so that upstream and downstream WSI grid boundaries 
coincide with global coarse grid boundaries. 

In the y and z directions, the WSI grid is evenly spaced near the store, and it 
is stretched near its outer boundary to provide reasonable overlap with the global 
coarse grid. Spacing is based on the parameter BAXIS which is defined to be the 
greater of 3.0 times the maximum store radius, excluding fins, or 1.5 times the 
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maximum store radius, including fins. In addition, a scaling factor is applied to 
BAXIS to limit its maximum value to 0.99 times the distance from wing plane to store 
centerline. Initially, 19 WSI grid points are evenly spaced in the y and z 
directions, extending a distance BAXIS inboard of, outboard of, above, and below 
the store centerline: 

y = B Q + Bjti (23) 

with 0.0 £ n £ 1.0, and 

z = C Q + CjC (24) 

with 0.0 £ £ £ 1.0. 

WSI grid spacing in the y and z directions is then modified to provide a 1-1/2 
global coarse grid mesh cell overlap between the WSI grid outer boundary and the 
global coarse grid inner boundary. The first three and last three WSI grid mesh 
cells in both the y and z directions are redefined using Eqs. 20 and 21. Polynomial 
coefficients are determined by requiring continuous first and second order grid 
metrics where original and redefined WSI grid regions abut, and by the overlap 
requirement. 

In some instances this procedure may generate a WSI grid which extends above 
(or below) the wing plane. When this occurs, polynomial coefficients are instead 
determined by requiring that the WSI grid outer boundary coincide with the wing 
plane itself. A new overlap region is then provided by extending the WSI grid an 
additional three mesh cells above (or below) the wing, again using Eq. 21. 

Polynomial coefficients for the additional mesh cells are determined by requiring 
continuous first and second order grid metrics at the wing plane, and by requiring a 
1-1/2 global coarse grid mesh cell overlap between the extended WSI grid boundary 
and the wing plane. 

The WSI grid inner boundary is defined by those y-z grid points which just 
enclose the store. This boundary is defined from just upstream of the store nose to 
the downstream grid boundary. 
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Coarse Cylindrical C-Grid 

Cylindrical and Cartesian coordinates are related by the following: 
y - y $ + r cos (9) 

(25) 

z = z $ + r sin (0). 

Angular spacing, 0, employs 24 evenly spaced mesh cells about each circumference. 
Coarse cylindrical C-grid physical coordinates x c and r c are nominally generated 
from computational coordinates £ and n via a series of conformal (Ref. 11) and 
shearing transformation: 

* “ ? ? MA X ^ 

n = n * n MAX 

x - - i 

Tlj = T1 

ri 2 = An + B tanh (Cn) 
n Loc = n 2 + ( mm ^ ^ n l " n 2^ 

n D = ^BODY + \0C ^ n 0UTER ^ " ^BODY (5)] 
n = tt - n D 

z = (x,r) = Jin (1 - cosh Z) where Z = (£,n) 

X c = X S IN G + "T" (x " An 2) 
r _ BAX I SO - 
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BAXISO is defined to be the greater of 4.0 times the store maximum radius, excluding 
fins, or the value of BAXIS before scaling. The scaling factor applied to BAXIS is 
then applied to BAXISO as well. Other parameters are L=Xj AIL -x SING and 
A =* 2it - Jtn(2). The variables n B0DY (?) and n 0UTER (?) shear the grid inner and 
outer boundaries to the specific store body and outer boundary shapes. The outer 
boundary is located a distance BAXIS from the store centerline and one quarter store 
body length upstream of the store nose. Forward of the store nose location it is 
ellipsoid in shape. The coefficients A, B, and C provide a clustering of grid 
points towards the store body surface. The parameters and are chosen so 
that 0.0 £ £ £ 1.0 along the store body and 0.0 £ n S 1.0 between the store body and 
the outer boundary. 

Singularity location x SING for the conformal mapping is internally 
calculated. For blunt noses it is located at the center of nose curvature. For 
sharp noses an iterative scheme was developed which determines x^q by requiring 
that d(n B Q D y)/d(5)=0 as £ goes to zero ( i . e . , at the nose). An option is available 
to override the internal calculation of x^q and to use a specified input value. 

All cases presented in this report use the default internal calculation. 


Coarse C-grid setup uses 40 % mesh cells between store nose and tail. A total 
of 5 n mesh cells is used between the store body and the grid outer boundary. A 
transition region is then defined which extends from the store body tail, £-p, to a 
location approximately one quarter store body length downstream of the tail, At 
the location a decoupled spacing Xg(£) and r G (n) is derived from the nominal 
spacing x c (5,n) and r c (5,q): 

dx c 

X D (?) = X C (? D ,0) + (? - ?) D ^ (? D ,0) (27) 

= r c (? D ,n) 


The coarse C-grid is required to transition from the nominal C-grid transformation 
x c and r c , at to the decoupled spacing x D and r G , at 


xU,q) = (1 - 3t 2 + 2t 3 ) x c (5,n) + (3t 2 - 2t 3 ) x D (?) 
r(l-,n) = (1 - 3t 2 + 2t 3 ) r c (£,n) + (3t 2 - 2t 3 ) r Q (n) 


(28) 
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where 


t . s ~ ;t 
5 d - 5 t 

Approximately 8 £ mesh cells are used in the transition region between and £q 
For subsonic freestreams, an additional 10 5 mesh cells and Eq. 14 are then used to 
stretch the coarse C-grid to downstream infinity, with next-to-last grid points in 
the global coarse grid and coarse C-grid being coincident. For supersonic 
freestreams, Eq. 22 is used to stretch the coarse C-grid to the global coarse grid 
downstream boundary, using that additional number (10 at most) of % mesh cells 
producing the smallest positive value of the constant A 3 . A maximum total of 60 
coarse C-grid \ mesh cells is permitted. 


Finally, the coarse C-grid is modified so that constant n coordinate lines 
closely resemble fin tip vortex streamlines. The streamline locations are derived 
(approximately) from the body thickness distribution and slender body theory: 

r t - R £ < 29 > 


which is integrated to give: 


r 


2,.. , . n 2 / „ \ „2,„ > 1 1/2 ^ 


[r‘(x 


FIN TIP LE 


) + R (x) - R (x 


FIN TIP LE 


) 3 ; x<x 


FIN TIP LE 


r(x) - ^ b FIN /2; x piN TIp LE * x S x piN Tip TE y30) 

^ 2(x fin tip te } + r2(x) " r2(x fin tip te^ 1/2; X>X FIN TIP TE 


V 




For arbitrary store body shape R(x), this fin tip vortex streamwise location will 
not pass smoothly along the fin tip itself. After eliminating the discontinuities 
in dr/dx at the fin tip leading and trailing edges, the desired grid line shape r(x) 
is obtained by introducing small bumps in the n transformation, along each £ = 
constant grid line, using Eq. 18. 
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Fine Cylindrical C-Grid 

The fine cylindrical C-grid uses the same angular spacing as the coarse 
cylindrical C-grid. In the | and n directions, mesh cell density in the fine C-grid 
is twice that of the coarse C-grid: 80 £ mesh cells from nose to tail, 
approximately 16 £ mesh cells from ^ to Sp and 10 n mesh cells between the store 
body and the grid outer boundary. A maximum total of 100 fine C-grid £ mesh cells 
is permitted. Rather than stretching the fine C-grid downstream boundary with Eqs. 
14 or 22, it is simply truncated one quarter store body length downstream of the 
tail (at the end of the transition region, £p). 

Wing Fine Grid 

As in the basic NASA/Grumman Transonic Wing-Body Code, wing fine grid arrays 
are set up at each global coarse grid wing spanwise station, and at two additional 
planes beyond the wing tip. For treatment of stores the outer boundary is redefined 
so as to exclude wing fine grid y-z points passing within the WSI grid inner 
boundary. This outer boundary modification is defined for the total upstream and 
downstream extent of the wing fine grid. 
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APPENDIX B 

FINITE DIFFERENCE APPROXIMATIONS 


Flow Equation Algorithm 

At each grid point, finite difference approximations are substituted for terms 

appearing in the governing flow equations, Eqs. 1 and 2. An upwind rotated 

difference scheme for the full potential equation (Ref. 12) extended for transonic 
small disturbance applications (Ref. 3) is used. Differencing in Cartesian 
coordinates x-y and z is analogous to differencing in cylindrical coordinates x-r 
and 0, respectively, and is omitted for brevity. 

Difference approximations in the physical domain are related to those in the 
computational space as follows: 

^x = + Vx 

<J> r = + Vr 

^xx = ♦gg^x + 2 ^ri^x n x + ♦nA + ^xx + ♦Ax \ 

^xr = ♦gg^r + ^gn ^x’V + ^r 1 ^ + ♦rm T 1 x 1 V + ♦g^xr + ♦Acr 

^rr = ^gg^r + ^gr^A + ♦nn T V + ^g^rr + ♦n T1 rT 

J 

Grid metrics are calculated with second order accurate central difference formulae, 
except at grid boundaries, where appropriate one-sided formulae are used. 

The parameter U^-4TV (where T, U, and V are from Eq. 2) determines whether a 
particular point is treated as subsonic or supersonic. At subsonic points, U^- 
4TV<0, centered finite difference expressions are given by: 
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*5 = 


* 1 - 1 . j 

, k ) /2A ? 



V 


♦u-i 

>k )/2An 



*0 = 


*i »j »k 

_ 1 )/2A0 




^i+l,j,k 

- — 

0) M , 

j ,k " (2 "u } 

♦i.j.k + * 1-1 


*5n 


k<l. 


j+l ,k**i+l ,j*l 

>k )/4454n 

6 

r nn 

= ^l.j+l.k 

-1*1. 

j, k -( 2 -i)*i 

, j ,k + *1 . j-1 . 

k )/{An) Z 

CD 

CD 

-e- 

= ,J,k+l 

- **i.i 

,k + *i ,j ,k 

_ 1 )/(A0) 2 



( 32 ) 


Note that terms with the + superscript indicate new potential values, while those 

o 

without it denote values from a previous grid sweep. At supersonic points, U - 
4TV > 0, the governing flow equation is recast into characteristic coordinates: 


(a 2 - u 2 - v 2 ) <t> .+ a 2 d> + - 6 + <b QQ = 0 (33) 

1 ’ y ss y nn r y r r 2 T 00 v 

where 

u 2 = \ [- (T-V) + /(T-V) 2 + U 2 ] 

v 2 = \ [+ (T-V) + /( T " V ) 2 + y2 ] 

a 2 = \ [+ (T+V) + J [T V) 2 + U 2 ] 

and 

♦ss ■ (“ 2 +xx * 2uv *xr * '' 2 *rr ) / (“ 2 * ^ 

*nn ■ < v2 »xx - 2uv *xr * "Xr)'^ * ^ 


50 



Contributions to the 4> term are centrally differenced using Eq. 32. Contri- 
butions to the (J) term are upwind differenced using formulae appropriate for local 
flow and grid properties. For example, at points where u£ + v£ >0 and where 

A I 

un x + vn r >0, the following are used: 




p 5n 


r nn 


(24> i,j,k "*1,j,k " 2 *1-1,j,k + ♦l-2 > j,k) / ( A 5) 2 
^1-l,j-l,k"*1 f j-l,k - *1-l # j,k + 2<f> i,j,k "*1 > j i k) /A5AT1 

Xj-l.k* ♦lJ^.k) / M 2 


•\ 


J 


(34) 


A periodic tridiagonal matrix solver was developed to facilitate a ring 
relaxation scheme. The relaxation process starts with the n-surface adjacent to the 
body (and sting). Each constant £,q-ring is relaxed in succession moving from the 
upstream centerline to the downstream boundary. The process is repeated for each q- 
surface and proceeds outward until the grid outer boundary is reached. 


The governing flow equation expressed in cylindrical coordinates is 
indeterminate along the grid centerline (r=0). Here, appropriate Cartesian 
formulations are employed. 


Far Field and Symmetry Plane Boundary Conditions 

For subsonic freestreams, various grid boundaries are stretched to correspond 
to infinity. The far field boundary condition imposed at x=-oo, y=±®, and z=±® is: 

4-0 ( 35 ) 

The far field boundary condition imposed at x=+® is: 


d) -d) = d) =0 

Y x Y xx Y xy 

(Cartesian) 

(36) 

A =6 = (b =0 

Y x y xx Y xr 

(cylindrical ) 
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For supersonic freestreams, grid boundaries are located at a finite distance 
from the configuration. At the upstream boundary the inflow boundary condition is: 


4> =* 4> x = o 


(37) 


At top, bottom, and side boundaries a radiation-type boundary condition 
(Ref. 13) is imposed.. The value of <f> at the boundary is found from interior 

A 

values using an extrapolation along local flow equation characteristics. The 
extrapolation is performed in an x-z plane directly above and below the 
configuration, in an x-y plane directly to the sides of the configuration, and in 
appropriately canted planes in between. The potential at the boundary is then set 
using an upwind expression for <|> x : 




which can be rearranged to give: 


(38) 


*i,j,k = (4(f) i-l,j,k ' *1-2, j,k + ( 2A ^ X ) ^ x ) /3 


(39) 


The radiation-type boundary condition is applied repretiti vely while marching 
downstream in the x direction. At the downstream boundary the outflow boundary 
condition is: 


(*x)i-l/2,j,k " Wi-3/2,j,k 


(40) 


or: 


x _ x , ^x)j-1 + ^x)j-2 / . . \ 

*i,j,k *1-1 ,j ,k (C x ) i + '*i-l ,j,k *i-2,j,k' 


(41) 


When a wing or fuselage is modelled, the symmetry conditions at y=0 are given 


by: 


d> - A =0 
y y ^xy 


(42) 
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Embedded Grid Boundaries 

The coarse cylindrical C-grid outer boundary is used to illustrate the Neumann- 

type embedded grid boundary condition. The location of each coarse C-grid outer 

boundary point in the WSI grid is determined once and appropriate interpolation 

parameters saved. The first step in updating a coarse C-grid outer boundary point 

is to obtain the values of <J> , <j> , and <J> at the corresponding point in the WSI 

x y z 

grid. These are then transformed to the coarse C-grid computational space: 


‘ *x x n * Vn * *z z n 

Then, using the relationship: 

" h till » k ~ h »J ~ 1 » k 

1 2An 

a dummy point just outside the coarse C-grid outer boundary is set: 


(43) 


(44) 


* 


D 

i ,j+l,k 


♦i 




+ 2Ag<j> Tl 


(45) 


This is sufficient treatment for most embedded grid boundaries. At the C-grid 

outer boundaries (and at the wing fine grid upstream boundary) it is also necessary 

to supply information for upwind differencing. First, upwind values of <J> , <f> , 

xx xy 

and <j> and centrally differenced values of <j> , <J> , and <j> are obtained at the 

xz yy yz zz 

corresponding WSI grid point. These are then transformed to the coarse C-grid 
computational space: 


^nn ^XX X T1 + ^yy^n + ^zz z n + ^xy x r^Ti + ^xz x ri z ri 

* 2 Vn z n * *x x nn + *y>m * *z z nn 

Then, using the relationship: 

<L n = ^i ,j+2,k~ 2< frj ,j+l ,k + ^i „i,k 

UN o 

(Anr 


(46) 


(47) 
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an additional dummy point outside the coarse C-grid is set: 


, j+2,k = 2 *1,J+l,k ' +1 f j,k + (Ati) * 


nn 


( 48 ) 


In addition to the Neumann-type boundary condition described above, a 
Dirichlet-type embedded grid boundary condition is also employed. This is used to 
set potentials in wing fine grid arrays beyond the wing tip, and to set the 
potential at a single grid point at the front of the coarse and fine cylindrical 
C-grids (in which the potential would otherwise vary by an indeterminate constant). 


Wing and Pylon Surfaces 

Wing surface boundary conditions are imposed via the <j> zz term of the governing 
flow equation (Ref. 1), as illustrated for wing lower surface treatment in the wing 
fine grid system. Using the difference expression: 

<t> z - (*° +1 - ( t ) k _ 1 )/2Az (49) 

a dummy point above the wing plane is set: 

<f> k+ 1 = + 2 Az<J> z (50) 

where <J> Z is given by Eq. 3. Dummy point values are then used in the appropriate 
finite difference expressions; for example: 

4> zz = ( 2 *k-1 " 2 +k + 2 Azc(, z )/(Az) 2 (51) 

(note that i and j subscripts have been omitted for convenience). In the wing/pylon 
junction, a combination of first and second order accurate expressions are used: 

(52) 

and the dummy point value is given by: 


♦z " Z 


^k+1 ^ 


'k-1 


k+1 


2Az 


Az 


<f>° +1 = (<f» M + 2<() k + 4 Az <J> Z ) /3 


(53) 
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In the wake, the Kutta condition is enforced using the wing circulation r (the jump 
in potential at the trailing edge) to set the dummy point value: 

C ■ - r < 54 > 

Pylon surface boundary conditions are imposed in a similar manner. In the wing 
fine grid a Z-scheme (Ref. 14) is employed, as illustrated for pylon outboard 
surface treatment. Using the following difference expression: 


*y ' *5 5 y * *n n y 


where, for ^>0, 


*n " (♦i.J*! - *i,j-l )/2AT1 


(55) 


a triadiagonal system of equations for dummy point values inboard of the pylon plane 
results: 


<2A T T 2A^ = 

T m) *i,j+i " m, *i±i,j+i ±<J) y 


where <j> is given by Eq. 5 (note that k subscripts have been omitted for 
convenience). The dummy point values are then used in the appropriate finite 
difference expressions. In the wing/pylon junction, combined first and second order 
accurate expressions are used: 


\ m \ [*^±1,3+1 " *i,j+l + *i,j-l " *iTl,j-l^ /2A5 
± ^ <f> i±l ,j “ *1,j + *i,j-l ‘ *1 t1,j-1^ /2A5 ^ 

= \ t^i.j+l ' *i,j-l ) /2Ati + (♦i.j " +? .j-i ) /At, 3 


A 


( 57 ) 


J 
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and the resulting tridiagonal system for dummy point values is given by: 


^ T 3A? *i,j-1 " *1 t1J-1 = 

\ T ^n } *i f >l " ^i±l , j+1 (58) 

+ **n"*i,j *i±i ,j ] * 

In the wake, the pylon circulation is used to set the dummy point value: 

tr (59) 
Store Body and Fin Surfaces 

The store body surface boundary condition, Eq. 7, is implemented using Eq. 31 
and the following finite difference expressions: 


V ,k> /2i5 

*11 ' (~ 3 *i,j,k * 4 *i , j+1 >k " *i,j+2,k) /24n 


to arrive at a tridiagonal system of equations for store body surface potentials: 


*i,j,k " ^x n x + ^r n r^ ^i+l,j,k *1-1 ,j ,k^ /2A? 
+ K n x + Vr) (^IJ+l.k ' ♦l,j + 2,k) /2Al1 


+ U cD n x + W ' t 3 (Vx + Vr) /2AT1 ^ 


Special care is required at the store nose, where for sharp noses the potential is 
singular and/or multivalued. To avoid numerical difficulties, the boundary 
condition at the nose itself is replaced by the requirement that the 
streamwise [%) derivative of the boundary condition be zero at body surface points 
adjacent to the nose: 


t( u . + *xK + + ‘VKl = 0 
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A store nose stagnation point potential value is also computed for use in flow 
computations at the point just upstream of the nose: 

+ <t> x = 0 (63) 

or 


■l.l.k 


(4<|) 


l,2,k 


- *1,3, k * T? U J /3 


(64) 


In a store body/fin junction the potential can also be set based on the fin 
boundary condition, Eq. 9, as illustrated for a body/fin upper surface junction. 
The difference expression: 



9 



+ 4(f) 


i , j ,k+l 


-4, iiJ>kf2 )/2r46 


is used to arrive at: 


(65) 


<J> 


+ 

i.j.k 


[4(|> i,j,k + l " *1,j,k+2 


2rA0(<J)g/r)]/3 


( 66 ) 


The tridiagonal matrix for store body surface potentials is modified at body/fin 
junction points to use the average of the Eq. 61 and Eq. 66 right hand side 
expressions . 


Store surface potentials are updated during each C-grid sweep using an under- 
relaxation factor of 0.5. Once all body surface potentials have been set, dummy 
potentials inside the body surface (for use by the rotated difference scheme) are 
set by equating two expressions for tf)^: 

x _ ^i »j+l ,k ~ ^1 ,j-1 ,k _ ~ 3( h,j,k +4< h ,j+l ,k ~^i ,j+2 ,k (67) 

% " 2An ■ 2Aq 

which yields: 


^i , j-1 ,k " 3< N,j,k 3< h ,j+l ,k +( h,j+2,k 


( 68 ) 


Store fin surfaces receive a lifting-surface treatment analogous to that for 
wing and pylon surfaces, as illustrated for a fin upper surface. Using the 
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following difference expression: 


7 4 >q = (* k +l ^k-l )/2rA0 (69) 

a dummy point below the fin plane is set: 

*k-l = ^k+l - 2rA0 (V r) (70) 

where <f>Q/r is given by Eq. 9 (note that i and j subscripts have been omitted for 
convenience). In the wake, the fin circulation is used to set the dummy point 
value: 



APPENDIX C 

CALCULATION OF FORCE AND MOMENT COEFFICIENTS 


Computed surface pressure and estimated skin friction coefficients are 
integrated to yield load distributions and force and moment coefficients. For the 
wing and fuselage this integration is identical to that found in the basic Transonic 
Wing-Body Code (Ref. 1). For cases where the wing finite difference boundary layer 
calculation is not activated, estimated wing upper and lower surface section skin 
friction coefficients are obtained in the same manner as for the fuselage, using the 
Prandtl-Schl ichting formula (Ref. 15) corrected for compressibility effects: 

c f,ave = (1 + °* 028M J 0.455/[log(Re)] 2 * 58 (72) 


where the Reynolds number is based on local chord. 


In general, the pylon surface will not be represented in its entirety in any 
one grid system. It therefore becomes necessary to piece together pylon coefficient 
contributions from several grid systems. First, contributions from that portion of 
the pylon located in the fine C-grid are computed. Section coefficients for the 
almost streamwise, constant ^-coordinate grid lines are defined as follows: 


c„ = 


c = 


m 


-?/J TE < c p, u - c p,i )d5 

^LE 

' C P,u- C P, 1» ( f-°- 25) 


'N 


ds 


c . = - - f TE (q n + C n ,)ds 
d s J^ LE v p,u x ,u p ,1 x,V 


(73) 


2 pTE 
'f = s J £LE 


Cj- ds 
f,ave 


y 
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where 


s = ds 
^LE 


"N 


ds = (n l + n$) 1/2 /(? x v K r \) d? 


(73) 

Cont'd 


J 


Contributions to total pylon force and moment coeffients in aircraft body axes are 
defined as follows: 


AF - ACX i + ACY n j + ACZ n k 
P P P P 


" A REF ^ =0 ^LE ^ P ’ u ^ + C P’ ini)dA 


AM p = ACMXpi + ACMY p j + ACMZ z k 




WreF ^ ^ ,C P.u ' C p,l )j X ( * ' *REF> dA 


(74) 




where 


n = n x i + n j + n z k 

X - i REF - (X - x REF ) i + yj + (z - z w1ng) k 
dA . l/(? x n r - 5 r n x ) d|dn 


A separate skin friction value is also calculated: 


AC 


F.P 


A REF 


fn=l [ KjE c dA 
J n=0 J ? LE C f ,ave aA 


(75) 


Next, contributions from that portion of the pylon surface located in the wing fine 
grid (but not in the fine C-grid) are considered. Contributions are computed using 
expressions analogous to Eqs. 73-75. Similarly, contributions from portions of the 
pylon surface in the global coarse grid and WSI grid are considered, as required, 
until the entire pylon surface has been accounted for. 
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Once the total pylon force and moment coefficients in aircraft body axes have 
been calculated, coefficients in aircraft stability axes are given by the following 
transformation: 


F p ■ S.p 1 + C y.P J + C L,P k ■ F P ' F P X 


M„ - C. „i + C m i + C„ „k - - M„ x aj 

p l,p m,p J n,p p p J 


(76) 


The skin friction coefficient remains unchanged. Finally, combined left- and right- 
hand pylon surface contributions to total configuration force and moment 
coefficients are given by: 


AC D " 2C D,p 


AC, = 2C, _ 
L L,p 


AC = 2C 
m m,p 


ACr- - 2C c „ 
F F,p 




(77) 


J 


Store body loads are calculated as follows. Cross-section force coefficients, 
in grid system axes and based on local diameter, are defined as: 


1 #2ir _ 

V + V + c z k - - 2 lo C p n de 


A 


where 


) (78) 


n = 


n i + n cos0j + n sin0k 

A I I 


J 


Store body viscous crossflow effects (Refs. 16, 17) are then estimated. First, a 
crossflow Reynolds number is defined at each store body cross-section in the fine 
C-grid: 
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(79) 


Re c = pV c 2r/pi 


where 


\ - v . C'° + °s > 2 * B s 2 ] V2 

The viscous crossflow drag coefficient for an infinite cylinder placed normal to the 
flow, c. , is then calculated as a function of Re , via a curve fit to 

U L. 

experimental data (see Ref. 16, Fig. 10). Next, the store fineness ratio is 
defined: 


L/d = ( X TAIL - X N0 $e )/(2*RMAXS) 


(80) 


The ratio of crossflow drag for a finite cylinder to that for an infinite 
cylinder, n c , is then calculated as a function of L/d, again via a curve fit to 
experimental data (see Ref. 16, Fig. 11). The estimated viscous crossflow drag 
coefficient value, referenced to freestream flow conditions, is given by: 





(81) 


This is broken down into cross-section viscous side and normal force coefficients: 


o 

II 

-hi 

I COS0 

y,c 1 c 


c -If 1 

s i n0 

z,c 1 c 1 



\ (82) 

where 

0 C = tan -1 [ (a + a $ )/3 s ] 

which can be added to cross-section coefficients calculated with Eq. 78. All loads, 
forces, and moments are subsequently calculated both with and without estimated 
viscous crossflow effects. 
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To calculate total store body force and moment coefficients the following 
quantities, in grid system axes, are first computed: 


F b - cy ♦ c? b j * cy 


■M 


'REFS s N0SE 


TAIL 


2rf ds 


M k = CMXi + CMYj + CMZk 


1 


X 


J* lfUl - 2rf x (x 

a REFS l REFS 


TAIL 

NOSE 


REFS 


)i ds 


A 




(83) 


J 


where ds is as defined for Eq. 73. Store body forces and moments in store body axes 
are then defined as: 

K - c v + c v + c v 

= CX b i + (CY b cos<J> s + CZ b sin4> s ) j + (CZ b cos<f> s - CY^in^Jk 

) (84) 

M b = CMX b i + CMY b j + CMZ b k 

= CMX b i + (CMY b cos<j> s + CMZ b sin<|> s ) j + (CMZ b cos<}> s - CMY b sin<|> s )k 


where <j> s is the store roll angle. Forces and moments in store stability axes are 
defined as: - 


F 


b 





« ^b " ^b x + a s^ " 8 s k ^ 




(85) 


M 


b 


C. ui ^ C l j ^ C uk 
A,b m,b J n,b 

M b - M b x [(a + a s )j - 6 s k] 


J 
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A separate skin friction value is also calculated using Eq. 72, with a Reynolds 
number based on store body length, and the following: 


P — b »W6t 
F»b ^reps f ,ave 


( 86 ) 


Store fin section coefficients are calculated using Eq. 73. To calculate total 
store fin force and moment coefficients the following quantities, in grid system 
axes, are first computed: 


F f = CX f i + CY f j + CZ f k 

X 


— C iP J* K u% + C n 1"1 )dA 
'REFS ’Voot ^LE P ’ u u P>1 1 


M r = CMX^i + 


CMY f j + CMZ f k 


= A T fn tiP ^ TE ( C D u " C d l)(' sin0 J + cos0k ) 

h REFS l REFS root ^LE P,U p, ‘ 

x(x - x REF$ )dA 

n - n i + (n cos0 - n fl sin0)j + (n fl cos0 + n sin0)k 


where 


x - X REFS = (x - x R p F<; )i + rcosSj + rsin8 k 


'REFS 


*\ 


)(8 7) 


J 


and dA is as defined for Eq. 74. Fin forces and moments in store body axes are then 
defined as: 


F f = CX f i + CY f j + CZ f k 




CX f i + (CY f cos«j> s + CZ f sin<f> s )j + (CZ f cos4> s - CY f sin<|> s )k 


( 88 ) 


M f = CMX f i + CMY f j + CMZ f k 


= CMX f i + (CMY f cos<J> s + CMZ f sin<J> $ ) j + (CMZ f cos<}> s - CMY f sin4> s )k 


J 
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Fin forces and moments in store stability axes are defined as: 


F f = C Djf i + C Yjf j + C L>f k 

= F f - F f x [(a + o $ )j - B g k] 

Mf = C ijf i + C m>f j + C njf k 

= M f - M f x [(a + a $ )j - 8 $ k] 

A separate skin friction value is also calculated: 



'N 


) (89) 


J 


(90) 


Total store forces and moments are obtained by summing contributions from the 
body and fins: 


F. = CX i + CY j + CZ k - F. + S F, 


M $ = CMX $ i + CMY $ j + CMZ $ k = + Z H f 

F s ■ C D, 5 ; + C V,sJ : * C L,s k ‘ F b + ZF f 

H s ‘ S.s 1 * C m, S j + C n,s k ‘ H b + E R f 


Cr = Cr 1, + Z Cr r 
F,s F,b F,f 


'N 


) (91) 




Finally, combined left- and right-hand store contributions to total configuration 
force and moment coefficients are given by: 


Ar _ , REFS r 
° ‘ D - 5 


- "REFS r 

iC L ' 2 ?WT C L,S 




/ (92) 
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GENERAL COMPUTER CODE DESCRIPTION 

The TSCLP computer code is operational on both CRAY XMP and CDC VPS-32 comput- 
ers. After main program execution, a second auxiliary program is used to generate 

the plotted output. Disk unit numbers 87 and 88 are used to pass information from 

the main program to the auxiliary plotting program. 

On the CRAY XMP, the second program module generates DISSPLA plotted output. On 
the CDC VPS-32, the second module does not actually call any plotting subroutines. 
Instead, it places all plotting subroutine call data on disc unit number 89 for 

transfer to the CDC NOS environment. On the NOS side, a small program module then 

generates Langley CALCOMP plotted output. 

Program execution requires approximately 780Kjg words of memory. In-core stor- 
age of data reduces I/O time and minimizes use of temporary disk storage units. In 
addition to the disk unit numbers discussed above, disk unit number 1 is used to 
store a copy of the input data set, and disk unit number 8 is used for quick geometry 
problem diagnosis output. Total CRAY execution time for various options are shown 
below. 


Case 

Geometry verification 
Isolated store body 
Store with 8 fins 
Wing/fuselage/pylon 
Wi ng/fuselage/pylon/store 
Wing/fuselage (supersonic) 


Iteration Count 
(Crude/Intermediate/Fi ne) 
0 / 0/0 
0/200/200 
0 / 200/200 
200 / 0/200 
200 / 200/200 
200 / 0/200 


Time 

(CPU seconds) 
12 
186 
224 
281 
400 
325 


Typical CDC execution time is approximately 19 CRUs for each Cray CPU second. 
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INPUT DATA FORMAT 


Excluding literal cards, all input data cards are punched in seven-field ten-digit 
format (7F10.0) . 

Wing/fuselage/pylon geometry is input in the aircraft coordinate system. 

Store body and fin geometry may be input in a separate store coordinate system. 
Parameters on Card 2-S are used to translate, rotate, and scale store geometry as 
desired within the aircraft coordinate system. 


Card 

Card 

Variable 


Number 

Col umn 

Name 

Description 

Card 1-A 

1-60 

TITLE 

Configuration or run title. 

Card 2-A 

1-10 

AMACH 

Freestream Mach number. 


11-20 

A0A 

Aircraft or isolated store angle-of-attack 
(degrees) . 


21-30 

RE 

Reynolds number X 1.0E-06 based on REFL (based on 
REFLS for isolated store). 


31-40 

AXITC 

Number of coarse grid iterations (global coarse 
grid only). Use 0.0 for isolated store. 


41-50 

AXITM 

Number of intermediate grid iterations (global 
coarse grid, wing/store interaction grid, and 
coarse C-grid) . 


51-60 

AXITF 

Number of fine grid iterations (global coarse 
grid, wing/store interaction grid, coarse C-grid, 
wing fine grid, and fine C-grid). 


61-70 

VISMOD 

=1.0 No boundary layer effects. 

=2.0 Wing boundary layer effects computed at end 
of inviscid analysis. 

=3.0 Wing invi scid/boundary layer interactive 
sol ution . 
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Card 3-A 1-10 

WING 

=0.0 No wing input. 

=+1.0 Wing input follows. 

=-1.0 Wing input follows. Wing included in grid 
setup, but not in solution process. 

11-20 

BODY 

=0.0 No fuselage input. 

=+1.0 Fuselage input follows. 

=-1.0 Fuselage input follows. Fuselage included 
in grid setup, but not in solution process. 

21-30 

STOR 

=0.0 No store input. 

=+1.0 Store input follows. 

=-1.0 Store input follows. Store included in grid 
setup, but not in solution process. 

NOTE: Omit card 4- 

-A for WING = 

BODY = 0.0 

Card 4-A 1-10 

REFA 

Aircraft reference area. 

11-20 

REFL 

Aircraft reference length. 

21-30 

REFX 

Aircraft reference X-position (for force and 
moment calculations). 
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NOTE: Omit card set -W for WING=0.0 


Card 

Card 

Variable 


Number 

Column 

Name 

Description 

Card 1-W 

1-10 

ASECT 

Number of streamwise sections defining wing 
planform (2.0 < ASECT < 20.0). 


11-20 

ANIN 

Number of ordinates defining each wing section 
(ANIN < 60.0). 


21-30 

ANOSW 

=0.0 Sharp nose wing sections. 
=1.0 Blunt nose wing sections. 


31-40 

ZWING 

Wing Z-position. 


41-50 

WS 

Wing C p distribution plot scaling per inch (typi- 
cally 0.4 or 0.8). 


51-60 

PYLS 

=0.0 No store pylon input. 

=1.0 Store pylon input follows. 

NOTE: Card set 2-W through 5-W is 

repeated ASECT times (sections will be 

extrapolated towards aircraft centerline as required) 

Card 2-W 

1-10 

XPL 

Wing section leading edge X-value. 


11-20 

YP 

Wing section Y-value. 


21-30 

XPT 

Wing section trailing edge X-value. 


31-40 

TWIST 

Wing section twist angle (degrees, positive lead- 
i ng edge up) . 


41-50 

AKODE 

■ =0.0 Section ordinates identical to preceding 
section (omit cards 4-W and 5-W). 

=1.0 New section definition expected on cards 4-W 
and 5-W. 

Card 3-W 

1-70 

XINW 

Wing section x/c coordinates (card 3-W defined 
only for first wing section, ANIN values). 

Card 4-W 

1-70 

YINU 

Wing section upper surface y/c coordinates (ANIN 
values). 

Card 5~W 

1-70 

YINL 

Wing section lower surface y/c coordinates (ANIN 


val ues) . 
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NOTE: Omit card set -B for BODY=0.0 


Card 

Card 

Variable 


Number 

Col umn 

Name 

Description 

Card 1-B 

1-10 

BKOD 

=2.0 Axisymmetric fuselage definition. 

=3.0 Quick-Geometry model fuselage definition. 


11-20 

BN0SE 

Fuselage nose X-value. 


21-30 

BTAIL 

Fuselage tail X-value. 


31-40 

BNIN 

Number of axisymmetric fuselage coordinates to be 
input (BNIN < 60.0) . 

(BKOD=2 . 0 only) 


41-50 

ANOSB 

=0.0 Sharp nose fuselage. 
=1.0 Blunt nose fuselage. 
(BK0D=2.0 only) 


51-60 

ZBODY 

Fuselage axis Z-value. 
(BK0D=2.0 only) 


61-70 

BS 

Fuselage C p plot scaling per inch (typically 
0.08). 

NOTE: Omit cards 2- 

■B and 3-B for BKOD=3.0 

Card 2-B 

1-70 

XINB 

Axisymmetric fuselage X-coordi nates (BNIN values). 

Card 3-B 

1-70 

RIN 

Axisymmetric fuselage radii (BNIN values). 

NOTE: Omit 

cards 4- 

B through 13-B 

for BK0D=2 . 0 

Card 4-B 

1-70 

VTITLE 

Quick-Geometry model title. 

Card 5-B 

1-10 

ACSM 

Number of distinct cross-section models (ACSM card 
sets 6-B and 7-B will follow). 

Card 6-B 

1-10 

ADUM 

Running count of current cross-section model 
(1-ACSM). 
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11-20 


AARC 


Number of arcs in current cross-section model 
(AARC Cards 7-B will follow). 

21-60 CTITLE Title or descriptor of current cross-section 

model . 


Card 7-B 

1-8 

ARCNAM 

Arc or component name. 


11-14 

ASHAPE 

Arc or component shape. 


21-28 

PNTNAM(l) 

Control point name for beginning of this arc. 


31-38 

PNTNAM(2) 

Control point name for termination of this arc. 


41-48 

PNTNAM(3) 

Slope control point name for this arc, if 
required. 

Card 8-B 

1-10 

ANTCSM 

Number of cross-section models to define entire 
fuselage (ANTCSM cards 9-B will follow). 

Card 9-B 

1-10 

ADUM 

Running count of current cross-section model 
(1-ANTCSM). 


11-20 

AMODEL 

Index corresponding to already defined cross- 
section models (between 1 and ACSM). 


21-30 

XCSMS1 

Starting X-value for current cross-section model. 


31-40 

XCSMS2 

Ending X-value for current cross-section model. 

Card 10-B 

1-10 

BLINE 

Number of fuselage line models to be defined 
(BLINE card sets 11-B and 12-B follow). 


11-20 

ALIAS 

Number of fuselage line models to be aliased 
(ALIAS cards 13-B follow). 

NOTE: Card 

set 11-B 

and 12-B is 

repeated BLINE times 

Card 1 1— B 

1-10 

BLSEG 

Number of segments defining fuselage line model. 


11 

BYORZ 

The letter Y or Z indicates which data definition 
is to follow. 


12-19 

BNAME 

Fuselage line name to be defined. 

Card 12-B 

1-4 

SSHAPE 

Segment shape. 


11-20 

0(1) 

X-value for beginning of segment. 
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21-30 

D(2) 

Y or Z value at D( 1) . 

31-40 

D(3) 

X-value for termination of segment. 

41-50 

D(4) 

Y or Z value at D(3) . 

51-60 

D(5) 

X-value for segment slope control point 

61-70 

D(6) 

Y or Z value at D(5) . 


NOTE: Card 13-B is repeated ALIAS times 


11 • 

BY0RZ 

The letter Y or Z indicates which data definition 
i s to fol low. 

12-19 

BNAME 

Fuselage line name to be defined. 

21 

AYORZ 

The letter Y or Z indicates which definition is to 
be used for aliasing. 

22-29 

ANAME 

Fuselage line name to which BNAME is aliased. 



NOTE: Omit card set -P for WING=0.0 or PYLS=0.0 


Card 

Card 

Variable 


Number 

Column Name 

Description 

Card 1-P 

1-10 

PSEC 

Number of streamwise sections defining pylon 
planform (2.0 < PSEC < 10.0). 

PSEC < 0.0; pylon included in grid setup but not in 
solution process. 


11-20 

PNIN 

Number of ordinates defining each pylon section 
(PNIN < 61.0). 


21-30 

PN0SE 

=0.0 Sharp nose pylon sections. 
=1.0 Blunt nose pylon sections. 


31-40 

YPYLS 

Pylon Y-position. 


41-50 

BETAP 

Pylon yaw angle (degrees, positive leading edge 
outboard) . 


51-60 

PS 

Pylon C distribution plot scaling per inch. 
P 

NOTE: Card 

set 2- 

P through 4-P is 

repeated PSEC times. Input sections in order of 

increasing 

distance from wing plane. Sections will be extrapolated towards 

wing plane 

and store (if present) as required. 

Card 2-P 

1-10 

XPLP 

Pylon section leading edge X-value. 


11-20 

YPP 

Pylon section Z-value. 


21-30 

XPTP 

Pylon section trailing edge X-value. 


31-40 

PKODE 

=0.0 Section ordinates identical to preceding 
section (omit card 4-P). 

=1.0 New section definition expected on card 4-P. 

Card 3-P 

1-70 

XINP 

Pylon section x/c coordinates (card 3-P defined 
only for first pylon section, PNIN values). 

Card 4-P 

1-70 

YINP 

Pylon section upper/lower surface y/c coordinates 
(PNIN values). 
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NOTE: Omit card set -$ for STOR=0.0 


Card 

Card 

Variable 


Number 

Col umn 

Name 

Descri pti on 

Card 1-S 

1-60 

TITLES 

Store identifying title. 

Card 2-S 

1-10 

SNOSE 

Store nose X-value in aircraft coordinate system. 


11-20 

STAIL 

Store tail X-value in aircraft coordinate system. 


21-30 

YST0R 

Store axis Y-value in aircraft coordinate system. 


31-40 

ZST0R 

Store axis Z-value in aircraft coordinate system. 


41-50 

ALPAS 

Store pitch angle with respect to aircraft 
(degrees, positive nose up). 


51-60 

BETAS 

Store yaw angle with respect to aircraft (degrees, 
positive nose outboard). ' 


61-70 

ROLLS 

Store roll angle (degrees, positive counterclock- 
wise looking upstream). 

NOTE: Remaining geometry may be 

input in a separate store coordinate system. 

Card 3-S 

1-10 

SNIN 

Number of coordinates defining axi symmetric store 
body (SNIN < 61.0). 


11-20 

ANOSES 

=0.0 Sharp nosed store. 
=1.0 Blunt nosed store. 


21-30 

ASING 

=0.0 C-grid singularity location XSING calculated 
i nternal ly . 

=1.0 C-grid singularity location XSING to be 
i nput. 


31-40 

XSING 

C-grid singularity location (ASING=1.0 only). 


41-50 

FINS 

=0.0 No store fins. 

=1.0 One set of store fins to be input. 
=2.0 Two sets of store fins to be input. 


51-60 

SS 

Store C p plot scaling per inch. 
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Card 4-S 

1-10 

REFAS 

Store reference area. 


11-20 

REFLS 

Store reference length. 


21-30 

REFXS 

Store reference X-position. 


31-40 

OMEGAS 

Store roll rate, p/V (rad/unit length). Positi 

oo 

counterclockwise looking upstream. 

Card 5-S 

1-70 

XINS 

Store body X-coordi nates (SNIN values). 

Card 6-S 

1-70 

RINS 

Store body radii (SNIN values). 



NOTE: Omit card set -F for STOR=0.0 or FINS=0.0 


Card 

Card 

Variable 


Number 

Column 

Name 

Descri pti on 

Card 1-F 

1-10 

FANG 

Number of fin angular positions (1.0 < FANG < 4.0). 


11-20 

FS 

Fin C distribution plot scaling per inch. 


21-30 

DFIN 

r 

=0.0 No fin deflections (omit cards 3-F). 
=1.0 Fins deflected (include cards 3-F). 

Card 2-F 

1-40 

ANGF 

Fin angular positions (degrees, positive counter- 
clockwise from outboard horizontal looking 
upstream). FANG values. Fins cannot be in same 
plane as pyl on . 


NOTE: Card set 3-F through 7-F is repeated FINS times (fore-fin input followed by 
aft-fin input for FINS=2.0) 

NOTE: Omit card 3-F for DFIN=0.0 

Card 3-F 1-40 DELTAF Fin deflection angles (degrees, positive leading 

edge counterclockwise looking upstream). FANG 
values . 

Card 4-F 1-10 FSEC Number of streamwi se sections defining fin plan- 

form (2.0 < FSEC < 10.0). 

FSEC < 0.0; fin included in grid setup but not in 
solution process. 

11-20 FNIN Number of ordinates defining each fin section 

( FNIN < 61.0). 

21-30 FN0SE =0.0 Sharp nose fin sections. 

=1.0 Blunt nose fin sections. 
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NOTE: Card set 5-F through 7-F is repeated FSEC times (sections will be extrapolated 
towards store axis as required) 

Card 5-F 1-10 XPLF Fin section leading edge X-value. 

11-20 YPF Fin section radial location. 

21-30 XPTF Fin section trailing edge X-value. 

31-40 FKODE =0.0 Section ordinates identical to preceding 

section (omit card 7-F). 

=1.0 New section definition expected on card 7-F. 

Card 6-F 1-70 XINF Fin section x/c coordinates (card 6-F defined only 

for first fin section, FNIN values). 

Card 7-F 1-70 YINF Fin section upper/lower surface y/c coordinates 

(FNIN values). 
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SAMPLE INPUT DATA SETS 


NTF 5 DEG TEST CONE 


0.60000 
0.00000 
ANALYSIS : 

0.00000 
0.00000 
BODY ALONE 

0.87489 
1 .00000 

0.00000 

200.00000 

200.00000 

1 .OOOOO 

0.00000 

100.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

15.00000 

60.11638 

0.00000 

8.74886 

0.00000 

50.00000 

0.00000 

0.00000 

0.04000 

0.00000 

0.00000 

5.00000 

10.00000 

20.00000 

30.00000 

40.00000 

45.00000 

48.00000 

100.00000 

52.00000 

55.00000 

60.00000 

70.00000 

80.00000 

90.00000 

0.00000 

0.43744 

0.87489 

1.74977 

2.62466 

3.49955 

3.93699 

4. 19946 
4.37443 

4.37443 

4.37443 

4.37443 

4.37443 

4.37443 

4.37443 


PATHFINDER 

I NOSECONE 






0.84000 

0.00000 

1.43750 

0.00000 

200.00000 

200.00000 

1 .OOOOO 

0.00000 

0.00000 

1.00000 





ANALYSIS : 

BODY ALONE 






0.00000 

100.00000 

0.00000 

0.00000 

6.00000 

0.00000 

0.00000 

41 .00000 

1 .00000 

0.00000 

0.00000 

0.00000 

0.20000 

0.00000 

162.29517 

14.37500 

50.00000 





0.00000 

0. 12500 

0.25000 

0.37500 

0.50000 

0.62500 

0.75000 

0.87500 

1 . OOOOO 

1 . 12500 

1.25000 

1.37500 

1.50000 

1.62500 

2.50000 

3.75000 

5.00000 

6.25000 

7.50000 

8.75000 

10.00000 

1 1 . 25000 

12.50000 

13.75000 

15.00000 

16.25000 

17.50000 

18.75000 

20.00000 

21.25000 

22.50000 

23 . 75000 

25.00000 

26.25000 

27.50000 

37 . 50000 

50.00000 

62.50000 

75.00000 

87 . 50000 

100.00000 


0.00000 

0.61250 

0.84750 

1.01500 

1 . 14620 

1.25000 

1.33620 

1.42000 

1 .50250 

1 .58370 

1.66500 

1 .74500 

1.82500 

1 .90370 

2.43370 

3. 13250 

3.76630 

4.33750 

4.85000 

5.30750 

5.71000 

6.06250 

6.36500 

6.61870 

6.82370 

6.98370 

7.09750 

7. 16500 

7 . 18750 

7. 18750 

7. 18750 

7. 18750 

7. 18750 

7. 18750 

7. 18750 

7. 18750 

7.18750 

7. 18750 

7.18750 

7. 18750 

7. 18750 



NACA RM L53F07 100 INCH BODY 


0.99000 

8.00000 

3. 167000 

0.00000 

0.00000 

1 .00000 

ANALYSIS : 

BODY ALONE 


0.00000 

100.00000 

0.00000 

30.00000 

1 . OOOOO 

0.00000 

78.53982 

10.00000 

50.00000 

0.00000 

0.00248 

0.00995 

0.34317 

0.60000 

0.90000 

12.00000 

18.00000 

24.00000 

54.00000 

60.00000 

66.00000 

96.00000 

100.00000 


0.00000 

0.01726 

0.03458 

0.20454 

0.27720 

0.35760 

2.36500 

3. 11200 

3 . 70800 

4.97200 

5.00000 

4.95600 

3.03100 

2.50000 



0.00000 

200.00000 

200.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.08000 

0.00000 

0.00000 

0.03939 
1 .50000 

30 . 00000 

72.00000 

0.08800 
3 . OOOOO 

36.00000 

78.00000 

0. 15504 
6 . OOOOO 

42.00000 

84 . 00000 

0.24042 

9.00000 

48.00000 

90.00000 

0.06883 
0.51360 
4. 15800 
4.82900 

0. 10292 
0.86640 
4.48900 
4.61000 

0. 13668 
1 .44600 
4.72000 
4.27400 

0. 17049 
1 .93600 
4.87600 
3.75400 



1 . ooooo 


GBU-15 CWW STORE WINGS CANARDS 


0.95000 6.00000 0.93750 0.00000 200.00000 200.00000 

0 . OOOOO 0 . OOOOO + 1 . OOOOO 

WINGS ON : CANARDS ON 


0.00000 

38.62500 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

51 .00000 

1 . OOOOO 

0.00000 

0.00000 

2.00000 

0.40000 

0.00000 

16.92000 

4 . 50000 

21.75000 





0.0 

0.03810 

0. 15230 

0.34210 

0.60675 

0.94520 

1.35620 

1 .83805 

2.38885 

3 . 00640 

3.68835 

4.43195 

5.23430 

6.09220 

7.00225 

7 . 96090 

8.96435 

10.00865 

11 .08965 

12.20310 

13.34460 

14.50970 

15.69370 

16.89200 

18.09985 

19.31250 

20.52515 

21.73300 

22.93130 

24. 11530 

25.28040 

26.42190 

27.53535 

28.61635 

29.66065 

30.66410 

31 .62275 

32.53280 

33.39070 

34. 19305 

34.93665 

35.61860 

36.23615 

36.78695 

37.26880 

37.67980 

38.01825 

38.28290 

38.47270 

38.58690 

38.62500 






0.0 

0.30505 

0.59540 

0.85520 

1 .07330 

1.29165 

1 .50395 

1 .68970 

1.82880 

1.87500 

1.87500 

1 .87500 

1.87500 

1 .87500 

1 . 87500 

1 .87500 

1 .87500 

1.87500 

1 .87500 

1.87500 

1 .88135 

1.98985 

2.08900 

2. 17145 

2.22805 

2.25000 

2.25000 

2 . 25000 

2.25000 

2 . 25000 

2.25000 

2.23365 

2. 19975 

2. 15050 

2.09390 

2.02345 

2 . OOOOO 

2 . OOOOO 

2.00000 

2 . OOOOO 

2 . OOOOO 

2 . OOOOO 

2 . OOOOO 

2 . OOOOO 

1 .99920 

1 .91185 

1 .83990 

1 .78365 

1 .74330 

1.71900 

1.71090 






4.00000 

0.40000 






45.00000 

135.00000 

225.00000 

315.00000 




+2.00000 

36.00000 

1 .00000 





2.72500 

1 .87500 

13.24000 

1 .00000 




0 . OOOOO 

0.0002b 

0.00111 

0.00247 

0.00432 

0.00660 

0.00924 

0.01215 

0.01527 

0.01847 

0 . 02000 

0 . 05000 

0. 10000 

0. 15000 

0 . 20000 

0.25000 

0 . 30000 

0.35000 

0 . 40000 

0.45000 

0.50000 

0.55000 

0 . 60000 

0.65000 

0.67857 

0 . 70000 

0.73000 

0.76000 

0.79000 

0.82000 

0.85000 

0 . 88000 

0.91000 

0.94000 

0.97000 

1 .00000 







0.00000 

0.00074 

0.00146 

0.00217 

0.00284 

0.00349 

0.00410 

0.00467 

0.00520 

0.00569 

0.00590 

0.00889 

0.01171 

0.01336 

0.01434 

0.01485 

0.01501 

0.01487 

0.01451 

0.01395 

0.01324 

0.01238 

0.01141 

0.01033 

0.00968 

0.00916 

0.00842 

0.00764 

0.00684 

0 . 00600 

0.00514 

0 . 00424 

0.00331 

0.00235 

0.00135 

0.00032 







3.72700 

3.61000 

4.53900 

1 .00000 




0 . OOOOO 

0.00320 

0.00632 

0.00924 

0.01187 

0.01415 

0.01600 

0.01736 

0.01819 

0,01847 

0.01847 

0.01847 

0.01847 

0.01847 

0.01847 

0.01847 

0.01847 

0.01847 

0.01847 

0.01847 

0.01847 

0.01847 

0.01847 

0.01847 

0.01847 

0.01839 

0.01800 

0.01729 

0.01626 

0.01488 

0.01321 

0.01 122 

0.00889 

0.00626 

0.00329 

0 . OOOOO 







+5.00000 

58.00000 

1 . OOOOO 





21 .41972 

2.00000 

38.46027 

1 . ooooo 




0.0 

0 . 00500 

0.01000 

0.01500 

0.02000 

0.02500 

0.03000 

0.03500 

0.04000 

0.04500 

0.05000 

0. 10000 

0 . 1 5000 

0 . 20000 

0.25000 

0 . 30000 

0.35000 

0 . 40000 

0.45000 

0 . 50000 

0.55000 

0.60000 

0.65000 

0 . 66000 

0.67000 

0 . 68000 

0.69000 

0 . 70000 

0.71000 

0.72000 

0.73000 

0.74000 

0.75000 

0.76000 

0.77000 

0.78000 

0.79000 

0 . 80000 

0.81000 

0.82000 

0.83000 

0.84000 

0.85000 

0 . 86000 

0.87000 

0.88000 

0.89000 

0 . 90000 

0.91000 

0.92000 

0.93000 

0.94000 

0.95000 

0 . 96000 
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678.418 

0.0 

1.0 



0.00048 

0.00644 

0.00972 

0.01330 

0.01832 

0.02499 

0.02969 

0.03340 

0.03647 

0.03909 

0.04132 

0.04320 

0.04476 

0.04600 

0.04694 

0.04759 

0.04799 

0.04814 

0.04808 

0.04781 

0.04737 

0.04676 

0.04601 

0.04512 

0.04412 

0.04302 

0.04182 

0.04054 

0.03920 

0.03780 

0.03490 

0.03193 

0.02895 

0.02597 

0.02446 

0.02295 

0.02141 

0.01985 

0.01825 

0.01662 

0.01496 

0.01327 

0.01 157 
- 0.00036 

0.00987 

0.00816 

0.00645 

0.00475 

0.00305 

0.00134 

0.00048 

- 0.00445 

- 0.00707 

- 0.00973 

- 0.01319 

- 0.01743 

- 0.02013 

- 0.02210 

- 0.02361 

- 0.02478 

- 0.02569 

- 0.02638 

- 0.02688 

- 0.02723 

- 0.02745 

- 0.02756 

- 0.02757 

- 0.02750 

- 0.02735 

- 0.02713 

- 0.02684 

- 0.02648 

- 0.02603 

- 0.02549 

- 0.02486 

- 0.02413 

- 0.02330 

- 0.02238 

- 0.02137 

- 0.02029 

- 0.01800 

- 0.01569 

- 0.01352 

- 0.01163 

- 0.01080 

- 0.01006 

- 0.00938 

- 0.00876 

- 0.00817 

- 0.00760 

- 0.00703 

- 0.00644 

- 0.00583 

- 0.00061 

- 0.00518 

- 0.00450 

- 0.00378 

- 0.00301 

- 0.00222 

- 0.00141 

502 . 209 

127.220 

712.723 

0.0 

1 .0 



0.00270 

0.00855 

0.01162 

0.01513 

0 . 02003 

0.02682 

0.03151 

0.03435 

0.03572 

0.03615 

0.03606 

0.03575 

0.03539 

0.03507 

0.03483 

0.03468 

0.03459 

0.03456 

0.03455 

0.03455 

0.03455 

0.03452 

0.03442 

0.03428 

0.03406 

0.03376 

0.03337 

0.03288 

0.03229 

0.03160 

0.02995 

0.02793 

0.02560 

0.02298 

0.02158 

0.02013 

0.01863 

0.01709 

0.01552 

0.01391 

0.01228 

0.01064 

0.00899 

- 0.00256 

0.00734 

0.00569 

0,00404 

0.00239 

0.00074 

- 0 . 0009 1 


160. 55203 


160.31653 


160.21167 


160. 16689 


85 



0.00270 

- 0.00238 

- 0.01860 

- 0.01912 

- 0.01564 

- 0.01531 

- 0.01490 

- 0.01492 

- 0.01521 

- 0.01527 

- 0.01418 

- 0.01375 

- 0.00968 

- 0.00880 

- 0.00285 

502.098 

127.240 

0.00295 

0 . 00573 

0.02212 

0.02422 

0.03184 

0.03254 

0.03441 

0.03437 

0.03217 

0.03146 

0.01989 

0.01840 

0.00874 

0.00709 

- 0.00280 

0.00295 

- 0.00016 

- 0.01269 

- 0.01364 

- 0.01569 

- 0.01570 

- 0.01524 

- 0.01517 

- 0.01531 

- 0.01541 

- 0.01449 

- 0.01404 

- 0.00991 

- 0.00903 

- 0.00309 

551 .010 

147.000 

0.00349 

0.00662 

0.02180 

0.02371 

0.03053 

0.03109 

0.03206 

0.03191 

0.02954 

0.02890 

0.01851 

0.01710 

0.00781 

0.00621 

- 0.00332 

0.00349 

0 . 00058 

- 0.01043 

- 0.01 121 

- 0.01382 

- 0.01426 

- 0.01707 

- 0.01748 

- 0.01918 

- 0.01931 

- 0.01738 

- 0.01682 

- 0.01 183 

- 0.01077 

- 0.00365 

598.080 

166.000 

0.00253 

0.00557 

0.01910 

0 . 0207 1 

0.02685 

0.02746 

0.02918 

0.02917 

0.02761 

0.02711 

0.01822 

0.01693 

0.00819 

0.00667 

- 0.00234 

0.00253 

- 0.00040 

- 0.01235 

- 0.01353 

- 0.01774 

- 0.01821 

- 0.02067 

- 0.02101 

- 0.02222 

- 0.02225 

- 0.01914 

- 0.01845 

- 0.01245 

- 0.01116 

- 0.00273 

650. 151 

187.000 

- 0.00139 

0.00167 

0.01555 

0.01722 

0.02377 

0.02449 

0.02704 

0.02719 

0.02667 

0.02637 

0.01972 

0.01866 

0.01094 

0.00954 

0.00163 

- 0.00139 

- 0.00444 

- 0.01598 

- 0.01707 

- 0.02085 

- 0.02129 

- 0.02341 

- 0.02365 

- 0.02396 

- 0.02381 


- 0.00521 

- 0.00783 

- 0.01897 

- 0.01838 

- 0.01509 

- 0.01497 

- 0.01495 

- 0.01499 

- 0.01536 

- 0.01535 

- 0.01324 

- 0.01266 

- 0.00787 

- 0.00691 

712.747 

0.0 

0.00745 

0.00946 

0.02602 

0.02756 

0.0331 1 

0.03357 

0.03425 

0.03403 

0.02976 

0.02772 

0.01686 

0.01527 

0.00544 

0.00379 

- 0.00213 

- 0.00413 

- 0.01435 

- 0.01486 

- 0.01566 

- 0.01560 

- 0.01512 

- 0.01511 

- 0.01557 

- 0.01564 

- 0.01352 

- 0.01292 

- 0.00811 

- 0.00715 

736.767 

0.0 

0.00828 

0.01020 

0.02535 

0.02676 

0.03151 

0.03181 

0.03169 

0.03140 

0.02740 

0.02561 

0.01563 

0.01412 

0.00461 

0.00302 

- 0.001 14 

- 0.00301 

- 0.01 180 

- 0.01226 

- 0.01473 

- 0.01521 

- 0.01786 

- 0.01820 

- 0.01941 

- 0.01929 

- 0.01618 

- 0.01546 

- 0.00965 

- 0.00849 

759.864 

0.0 

0 . 007 1 5 

0.00897 

0.02209 

0.02329 

0.02798 

0.02840 

0.02908 

0.02892 

0.02592 

0.02444 

0.01557 

0.01416 

0.00516 

0.00366 

- 0.00221 

- 0.00410 

- 0.01453 

- 0.01537 

- 0.01866 

- 0.01909 

- 0.02131 

- 0.02158 

- 0.02213 

- 0.02177 

- 0.01767 

- 0.01681 

- 0.00982 

- 0.00843 

785.392 

0.0 

0.00328 

0.00511 

0.01866 

0.01993 

0.02512 

0.02566 

0.02727 

0.02728 

0.02557 

0.02452 

0.01752 

0.01632 

0.00814 

0.00678 

- 0.00621 

- 0.00802 

- 0.01797 

- 0.01872 

- 0.02171 

- 0.02210 

- 0.02383 

- 0.02397 

- 0.02334 

- 0.02263 


- 0.01 124 

- 0.01510 

- 0.01759 

- 0.01680 

- 0.01491 

- 0.01489 

- 0.01504 

- 0.01509 

- 0.01519 

- 0.01482 

- 0.01201 

- 0.01129 

- 0.00591 

- 0.00490 

1 .0 

0.01237 

0.01650 

0.02888 

0 . 03003 

0.03393 

0.03418 

0.03371 

0.03329 

0.02537 

0.02275 

0.01367 

0.01204 

0.00214 

0.00050 

- 0.00667 

- 0.00961 

- 0.01523 

- 0.01547 

- 0.01551 

- 0.01542 

- 0.01513 

- 0.01516 

- 0.01551 

- 0.01515 

- 0.01226 

- 0.01154 

- 0 .00615 

- 0.00514 

1.0 

0.01293 

0.01673 

0.02797 

0.02899 

0.03201 

0.03211 

0.03104 

0.03061 

0.02352 

0.02115 

0.01258 

0.01 100 

0.00144 

- 0.00015 

- 0.00529 

- 0.00785 

- 0.01265 

- 0.01302 

- 0.01570 

- 0.01617 

- 0.01851 

- 0.01878 

- 0.01892 

- 0.01829 

- 0.01466 

- 0.01379 

- 0.00729 

- 0.00608 

1 .0 

0.01 141 

0.01476 

0.02436 

0.02530 

0.02873 

0.02897 

0.02869 

0.02840 

0.02267 

0.02059 

0.01270 

0.01122 

0 . 002 1 6 

0 . 00066 

- 0.00637 

- 0.00912 

- 0.01608 

- 0.01670 

- 0.01951 

- 0.01992 

- 0.02182 

- 0.02200 

- 0.02116 

- 0.02029 

- 0.01586 

- 0.01481 

- 0.00702 

- 0.00558 

1 .0 

0.00765 

0.01 107 

0.02105 

0.02206 

0.02613 

0.02651 

0.02722 

0.02710 

0.02320 

0.02161 

0.01504 

0.01371 

0.00546 

0 . 004 1 7 

- 0.01027 

- 0.01292 

- 0.01934 

- 0.01989 

- 0.02248 

- 0.02282 

- 0.02405 

- 0.02408 

- 0.02168 

- 0.02045 


- 0.01730 

- 0.01614 

- 0.01488 

- 0.01515 

- 0.01453 

- 0.01051 

- 0.00387 

160.21803 

0.01960 
0.03100 
0.03434 
0.03278 
0.02135 
0.01040 
- 0.001 15 

- 0.01141 

- 0.01561 

- 0.01533 

- 0.01523 

- 0.01485 

- 0.01075 

- 0.00411 

159.90466 

0.01954 

0.02984 

0.03213 

0.03011 

0.01986 

0.00941 

- 0.00173 

- 0.00937 

- 0.01341 

- 0.01663 

- 0.01900 

- 0.01787 

- 0.01284 

- 0.00486 

159.33997 

0.01718 

0.02612 

0.02911 

0.02804 

0.01944 

0.00971 

- 0.00084 

- 0.01093 

- 0.01724 

- 0.02030 

- 0.02214 

- 0.01975 

- 0.01367 

- 0.00415 

158.37964 

0.01356 

0.02296 

0.02681 

0.02692 

0.02070 

0.01234 

0.00290 

- 0.01464 

- 0.02039 

- 0.02314 

- 0.02405 

- 0.01973 
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-0.01893 

-0.01805 

-0.01707 

-0.01598 

-0.01478 

-0.01347 

-0.01206 


-0.01055 

-0.00897 

-0.00733 

-0.00565 

-0.00395 

-0.00224 

-0.00053 


0.00116 

699.743 

-0.00675 

207.000 

-0.00368 

808.049 

-0.00195 

0.0 

-0.00006 

1.0 

0.00256 

0.00615 

0.00874 

157.43750 

0.01084 

0.01262 

0.01416 

0.01552 

0.01674 

0.01783 

0.01883 


0.01974 

0.02057 

0.02133 

0.02202 

0.02263 

0.02318 

0.02366 


0.02408 

0.02444 

0.02472 

0.02495 

0.02512 

0.02522 

0.02527 


0.02526 

0.02519 

0.02488 

0.02434 

0.02359 

0.02261 

0.02203 


0.02140 

0.02070 

0.01994 

0.01912 

0.01822 

0.01727 

0.01624 


0.01517 

0.01407 

0.01296 

0.01188 

0.01086 

0.00992 

0 . 00902 


0.00813 

-0.00675 

-0.00985 

-0.01167 

-0.01348 

-0.01565 

-0.01818 

-0.01979 


-0.02100 

-0.02195 

-0.02271 

-0.02333 

-0.02386 

-0.02433 

-0.02476 


-0.02516 

-0.02552 

-0.02584 

-0.02613 

-0.02637 

-0.02655 

-0.02669 


-0.02677 

-0.02679 

-0.02676 

-0.02667 

-0.02653 

-0.02634 

-0.02609 


-0.02579 

-0.02544 

-0.02459 

-0.02351 

-0.02218 

-0.02058 

-0.01965 


-0.01863 

-0.01751 

-0.01628 

-0.01492 

-0.01345 

-0.01185 

-0.01015 


-0.00837 

-0.00652 

-0.00461 

-0.00267 

-0.00069 

0.00132 

0.00335 


0.00537 

754.789 

-0.01701 

229.200 

-0.01392 

833. 199 
-0.01196 

0.0 

-0.00996 

1 .0 
-0.00718 

-0.00327 

-0 . 00049 

156.39172 

0.00182 

0 . 0038 1 

0.00554 

0.00708 

0.00849 

0.00973 

0.01092 


0.01202 

0.01307 

0.01407 

0.01505 

0.01593 

0.01680 

0.01763 


0.01841 

0.01918 

0.01984 

0.02049 

0.02110 

0.02162 

0.02211 


0.02256 

0.02293 

0.02356 

0.02400 

0.02434 

0.02452 

0.02458 


0.02462 

0.02461 

0.02457 

0.02448 

0.02431 

0.02409 

0.02371 


0.02327 

0.02274 

0.02219 

0.02164 

0.02120 

0.02093 

0.02074 


0.02057 

-0.01701 

-0.02021 

-0.02212 

-0.02393 

-0.02595 

-0.02825 

-0.02965 


-0.03061 

-0.03129 

-0.03178 

-0.03216 

-0.03251 

-0.03283 

-0.03313 


-0.03341 

-0.03362 

-0.03375 

-0.03385 

-0.03382 

-0.03369 

-0.03349 


-0.03320 

-0.03280 

-0.03237 

-0.03184 

-0.03128 

-0.03067 

-0.03000 


-0.02929 

-0.02856 

-0.02698 

-0.02519 

-0.02314 

*-0.02083 

-0.01950 


-0.01806 

-0.01648 

-0.01477 

-0.01289 

-0.01090 

-0.00875 

-0.00649 


-0.00420 

-0.00183 

0 . 00060 

0.00304 

0.00555 

0.00814 

0.01078 


0.01343 
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F-14 FUSELAGE MODEL 
8 . 

1. 2. NOSE 


BLO 

ELLI 

BCL 

MHB 

SCPLO 

BUP 

ELLI 

MHB 

TCL 

SCPUP 

2. 

4 . 

NOSE WITH 

1 CANOPY 


BLO 

ELLI 

BCL 

MHB 

SCPLO 

BSI 

LINE 

MHB 

FUP 


BUP 

ELLI 

FUP 

CREASE 

SCPUP 

CAN 

ELLI 

CREASE 

TCL 

SCPCAN 

3. 

7. 

INLET 



BLO 

ELLI 

BCL 

MHB 

SCPLO 

BSI 

LINE 

MHB 

FLO 


NACLO 

LINE 

FLO 

INLTLO 


NACSI 

LINE 

INLTLO 

INLTUP 


NACUP 

LINE 

INLTUP 

FUP 


BUP 

ELLI 

FUP 

CREASE 

SCPUP 

CAN 

ELLI 

CREASE 

TCL 

SCPCAN 

4. 

6. 

INLET TO 

NACELLE 

FOWARD FAIRING 

BLO 

ELLI 

BCL 

FLO 

SCPLO 

NACLO 

LINE 

FLO 

INLTLO 


NACSI 

LINE 

INLTLO 

INLTUP 


NACUP 

LINE 

INLTUP 

FUP 


BUP 

ELLI 

FUP 

CREASE 

SCPUP 

CAN 

ELLI 

CREASE 

TCL 

SCPCAN 

5. 

6. 

INLET TO 

NACELLE 

AFT FAIRING 

BLO 

LINE 

BCL 

FLO 


NACLO 

LINE 

FLO 

INLTLO 


NACSI 

LINE 

INLTLO 

INLTUP 


NACUP 

LINE 

INLTUP 

FUP 


BUP 

ELLI 

FUP 

CREASE 

SCPUP 

CAN 

ELLI 

CREASE 

TCL 

SCPCAN 

6. 

7. 

FOWARD NACELLE 


BLO 

LINE 

BCL 

FLO 


NACLO 

LINE 

FLO 

INLTLO 



87 



NACLS 

ELLI 

INLTLO 

MHB 

SCPLO 



NACSI 

LINE 

MHB 

INLTUP 




NACUP 

LINE 

INLTUP 

FUP 




BUP 

ELLI 

FUP 

CREASE 

SCPUP 



CAN 

ELLI 

CREASE 

TCL 

SCPCAN 



7. 

6. 

MID NACELLE 




BLO 

LINE 

BCL 

FLO 




NACLO 

LINE 

FLO 

INLTLO 




NACLS 

ELLI 

INLTLO 

MHB 

SCPLO 



NACSI 

LINE 

MHB 

INLTUP 




NACUP 

LINE 

INLTUP 

FUP 




BUP 

ELLI 

FUP 

TCL 

SCPUP 



8. 

6 . 

AFT NACELLE 




BLO 

LINE 

BCL 

FLO 




NACLO 

LINE 

FLO 

INLTLO 




NACLS 

ELLI 

INLTLO 

CREASE 

SCPCAN 



NACSI 

LINE 

CREASE 

INLTUP 




NACUP 

LINE 

INLTUP 

FUP 




BUP 

8. 

1 . 

LINE 

FUP 

TCL 




1 . 

93.0 

205.0 




2. 

2. 

205.0 

352.0 




3. 

3. 

352.0 

414.5 




4. 

4. 

414.5 

425.0 




5. 

5. 

425.0 

433.0 




6. 

6. 

433.0 

497.0 




7. 

7. 

497.0 

650.0 




8. 

8. 

650.0 

780.0 




22.0 

2.0 






10.0 

2TCL 






ELLX 

93.0 

131 .45 

205.0 

171 .0 

127.5 

152.0 

LINE 

205.0 

171.0 

241.0 

191 .65 



ELLY 

241 .0 

191.65 

313.0 

203.0 

274.0 

203.0 

ELLY 

313.0 

203.0 

352.0 

199.5 

331 .0 

203.0 

LINE 

352.0 

199.5 

497.0 

174.5 



LINE 

497.0 

174.5 

546.0 

169.5 



LINE 

546.0 

169.5 

600.0 

165.0 



LINE 

600.0 

165.0 

650.0 

160.5 



LINE 

650.0 

160.5 

697.0 

165.5 



LINE 

697.0 

165.5 

780.0 

162.75 



10.0 

ZBCL 






ELLX 

93.0 

131 .45 

176.0 

116.5 

129.5 

116.5 

LINE 

176.0 

116.5 

205.0 

118.0 



LINE 

205.0 

1 18.0 

352.0 

121.75 



LINE 

352.0 

121 .75 

425.0 

123.43 



LINE 

425.0 

123.43 

433.0 

117.0 



LINE 

433.0 

117.0 

497.0 

107.8 



LINE 

497.0 

107.8 

546.0 

104.7 



LINE 

546.0 

104.7 

600.0 

103.0 



LINE 

600.0 

103.0 

650.0 

104.0 



LINE 

650.0 

104.0 

780.0 

113.0 



7.0 

YMHB 






ELLY 

93.0 

0.0 

352.0 

31 .0 

164.0 

31 .0 

LINE 

352.0 

31 .0 

414.5 

29.5 



LINE 

414.5 

29.5 

433.0 

75.0 



LINE 

433.0 

75.0 

497.0 

74.0 



LINE 

497.0 

74.0 

546.0 

73.0 



LINE 

546.0 

73.0 

600.0 

76.5 



LINE 

600.0 

76.5 

650.0 

79.0 



9.0 

ZMHB 






LINE 

93.0 

131 .45 

164.0 

138.0 



LINE 

164.0 

138.0 

205.0 

138.5 



LINE 

205.0 

138.5 

352.0 

130.5 



LINE 

352.0 

130.5 

414.5 

127.93 



LINE 

414.5 

127.93 

433.0 

119.0 



LINE 

433.0 

1 19.0 

497.0 

125.5 



LINE 

497.0 

125.5 

546.0 

139.0 



LINE 

546.0 

139.0 

600.0 

139.5 



LINE 

600.0 

139.5 

650.0 

140.0 



9.0 

YCREASE 






LINE 

205.0 

0.0 

221 .0 

15.5 



LINE 

221 .0 

15.5 

283.0 

16.0 



LINE 

283.0 

16.0 

328.0 

16.5 




LINE 

328.0 

16.5 

LINE 

352.0 

16.0 

LINE 

405.0 

13.0 

LINE 

497.0 

0.0 

LINE 

650.0 

79.0 

LINE 

697.0 

80.5 

9.0 

ZCREASE 


LINE 

205.0 

171 .0 

LINE 

221 .0 

169.0 

LINE 

283.0 

179.0 

LINE 

328.0 

183.0 

LINE 

352.0 

183.5 

LINE 

405.0 

182.0 

LINE 

497.0 

174.5 

LINE 

650.0 

140.0 

LINE 

697.0 

141 .0 

7.0 

YFUP 


ELLY 

93.0 

0.0 

LINE 

352.0 

31 .0 

LINE 

497.0 

19.5 

LINE 

546.0 

14.0 

LINE 

600.0 

11.0 

LINE 

650.0 

0.0 

LINE 

697.0 

49.5 

8.0 

ZFUP 


LINE 

93.0 

131 .45 

LINE 

205.0 

138.5 

LINE 

352.0 

155.0 

LINE 

497.0 

163.0 

LINE 

546.0 

163.5 

LINE 

600.0 

163.0 

LINE 

650.0 

160.5 

LINE 

697.0 

165.5 

6.0 

YSCPUP 


ELLY 

93.0 

0.0 

LINE 

352.0 

31 .0 

LINE 

405.0 

23.5 

LINE 

497.0 

14.0 

LINE 

546.0 

7.5 

LINE 

600.0 

5.0 

7.0 

ZSCPUP 


ELLX 

93.0 

131.45 

LINE 

205.0 

171 .0 

LINE 

221 .0 

158.0 

LINE 

352.0 

173.0 

LINE 

497.0 

174.5 

LINE 

546.0 

169.5 

LINE 

600.0 

165.0 

9.0 

YSCPCAN 


LINE 

205.0 

0.0 

LINE 

221.0 

12.5 

LINE 

283.0 

16.0 

LINE 

328.0 

16.5 

LINE 

352.0 

16.0 

LINE 

405.0 

9.5 

LINE 

497.0 

0.0 

LINE 

650.0 

81 .5 

LINE 

697.0 

82.5 

7.0 

ZSCPCAN 


LINE 

205.0 

171.0 

ELLY 

241.0 

191.65 

ELLY 

313.0 

203.0 

LINE 

352.0 

199.5 

LINE 

497.0 

174.5 

LINE 

650.0 

105.5 

LINE 

697.0 

1 10.0 

9.0 

ZSCPLO 


ELLX 

93.0 

131.45 

LINE 

176.0 

116.5 

LINE 

205.0 

118.0 

LINE 

352.0 

121 .75 

LINE 

425.0 

123.43 

LINE 

433.0 

119.0 


352.0 16.0 

405.0 13.0 

497.0 0.0 

650.0 79.0 

697.0 80.5 

780.0 81.0 

221.0 169.0 

283.0 179.0 

328.0 183.0 

352.0 183.5 

405.0 182.0 

497.0 174.5 

650.0 140.0 

697.0 140.0 

780.0 140.0 

352.0 31 .0 164.0 

497.0 19.5 

546.0 14.0 

600.0 11.0 

650.0 0.0 

697.0 49.5 

780.0 51.0 

205.0 138.5 

352.0 155.0 

497.0 163.0 

546.0 163.5 

600.0 163.0 

650.0 160.5 

697.0 165.5 

780.0 162.75 

352.0 31.0 164.0 

405.0 23.5 

497.0 14.0 

546.0 7,5 

600.0 5.0 

650.0 0.0 

205.0 171.0 127.5 

221.0 158.0 

352.0 173.0 

497.0 174.5 

546.0 169.5 

600.0 165.0 

650.0 160.5 

221.0 12.5 

283.0 16.0 

328.0 16.5 

352.0 16.0 

405.0 9.5 

497.0 0.0 

650.0 81.5 

697.0 82.5 

780.0 83.0 

241.0 191.65 

313.0 203.0 274.0 

352.0 199.5 331.0 

497.0 174.5 

650.0 105.5 

697.0 110.0 

780.0 116.0 

176.0 116.5 129.5 

205.0 118.0 

352.0 121.75 

425.0 123.43 

433.0 119.0 

497.0 113.0 


31 .0 


31 .0 


152.0 


203.0 

203.0 


116.5 



LINE 

497.0 

113.0 

546.0 

110.0 

LINE 

546.0 

110.0 

600.0 

105.0 

LINE 

600.0 

105.0 

650.0 

105.5 

8.0 

YSCPLO 




ELLY 

93.0 

0.0 

352.0 

31 .0 

LINE 

352.0 

31 .0 

414.5 

29.5 

LINE 

414.5 

29.5 

425.0 

39.79 

LINE 

425.0 

39.79 

433.0 

75.0 

LINE 

433.0 

75.0 

497.0 

76.2 

LINE 

497.0 

76.2 

546.0 

77.7 

LINE 

546.0 

77.7 

600.0 

79.0 

LINE 

600.0 

79.0 

650.0 

81 .5 

1 .0 

YMAPAXIS 




LINE 

93.0 

0.0 

780.0 

0.0 

2.0 

ZMAPAXIS 




LINE 

93.0 

131.45 

352.0 

146.0 

LINE 

352.0 

146.0 

780.0 

147.0 

9.0 

YFLO 




LINE 

352.0 

31 .0 

414.5 

29.5 

LINE 

414.5 

29.5 

425.0 

39.79 

LINE 

425.0 

39.79 

433.0 

42.0 

LINE 

433.0 

42.0 

497.0 

46.3 

LINE 

497.0 

46.3 

546.0 

47.5 

LINE 

546.0 

47.5 

600.0 

48.5 

LINE 

600.0 

48.5 

650.0 

49.5 

LINE 

650.0 

49.5 

697.0 

49.75 

LINE 

697.0 

49.75 

780.0 

51 .0 

9.0 

ZFLO 




LINE 

352.0 

150.0 

414.5 

127.93 

LINE 

414.5 

127.93 

425.0 

123.43 

LINE 

425.0 

123.43 

433.0 

117.0 

LINE 

433.0 

1 17.0 

497.0 

107.8 

LINE 

497.0 

107.8 

546.0 

104.7 

LINE 

546.0 

104.7 

600.0 

103.0 

LINE 

60C.0 

103.0 

650.0 

104.0 

LINE 

650.0 

104.0 

697.0 

107.0 

LINE 

697.0 

107.0 

780.0 

113.0 

7.0 

YINLTUP 




LINE 

352.0 

66.5 

433.0 

66.5 

LINE 

433.0 

66.5 

497.0 

66.5 

LINE 

497.0 

66.5 

546.0 

68.0 

LINE 

546.0 

68.0 

600.0 

74.5 

LINE 

600.0 

74 . 5 

650.0 

77.0 

LINE 

650.0 

77.0 

697.0 

78.5 

LINE 

697.0 

78.5 

780.0 

80.0 

7.0 

ZINLTUP 




LINE 

352.0 

162.5 

433.0 

169.0 

LINE 

433.0 

169.0 

497.0 

171 .5 

LINE 

497.0 

171.5 

546.0 

170.0 

LINE 

546.0 

170.0 

600.0 

165.0 

LINE 

600.0 

165.0 

650.0 

162.0 

LINE 

650.0 

162.0 

697.0 

160.0 

LINE 

697.0 

160.0 

780.0 

155.5 

7.0 

YINLTLO 




LINE 

352.0 

67.0 

433.0 

75.0 

LINE 

433.0 

75.0 

497.0 

65.5 

LINE 

497.0 

66.5 

546.0 

61 .0 

LINE 

546.0 

61.0 

600.0 

57.0 

LINE 

600.0 

57.0 

650.0 

57.5 

LINE 

650.0 

57.5 

697.0 

58.0 

LINE 

697.0 

58.0 

780.0 

59.0 

7.0 

ZINLTLO 




LINE 

352.0 

156.5 

433.0 

1 19.0 

LINE 

433.0 

1 19.0 

497.0 

111.0 

LINE 

497.0 

111.0 

546.0 

107.0 

LINE 

546.0 

107.0 

600.0 

104.0 

LINE 

600.0 

104.0 

650.0 

104.0 

LINE 

650.0 

104.0 

697.0 

108.0 

LINE 

697.0 

108.0 

780.0 

114.0 


YTCL 

YMAPAXIS 




YBCL 

YMAPAXIS 
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200.0 


1.0 


SC 3 DEMONSTRATION WING - BASIC LEADING EDGE 


1 .620 

12.000 

2.4 

200 .0 0.0 

1 .O 

0.0 

0.0 


342. 1 1 

23.840 

16.701 


11.0 

30.0 

1 .0 

0.0 0.4 


0.0 

0.0 

23.840 

0.0 

1 .0 



0.0 

0.00147 

0.00586 

0.01317 

0.02338 

0.03645 

0.05235 

0.07102 

0.09242 

0. 11649 

0. 14314 

0 . 17231 

0.20391 

0.23784 

0.27400 

0.31230 

0.35261 

0.39483 

0.43881 

0.48445 

0.53159 

0.5801 1 

0.62986 

0.68070 

0.73247 

0.78503 

0.83822 

0.89188 

0.94586 

1 .00000 






0.0 

0.00187 

0.00367 

0 . 0054 1 

0 . 00709 

0.00870 

0.01023 

0.01169 

0.01306 

0.01434 

0.01552 

0.01659 

0.01754 

0.01836 

0.01903 

0.01954 

0.01987 

0.02000 

0.01991 

0.01959 

0.01900 

0.01814 

0.01699 

0.01554 

0.01378 

0.01171 

0.00933 

0.00665 

0.00367 

0.00040 






0.0 

- 0.00187 

- 0.00367 

- 0.00541 

- 0.00709 

- 0.00870 

- 0.01023 

- 0.01 169 

- 0.01306 

- 0.01434 

- 0.01552 

- 0.01659 

- 0.01754 

- 0.01836 

- 0.01903 

- 0.01954 

- 0.01987 

- 0.02000 

- 0.01991 

- 0.01959 

- 0.01900 

- 0.01814 

- 0.01699 

- 0.01554 

- 0.01378 

- 0.01171 

- 0.00933 

- 0.00665 

- 0.00367 

- 0.00040 






3. 151 

1 .470 

24. 135 

0.0 

1 .0 



- 0.00883 

- 0.00607 

- 0.00286 

0.00064 

0.00425 

0.00790 

0.01 152 

0.01503 

0.01839 

0.02153 

0.02436 

0.02691 

0.02913 

0.03104 

0.03263 

0.03386 

0.03476 

0.03534 

0.03557 

0.03547 

0.03502 

0.03424 

0.03313 

0.03167 

0.02989 

0.02778 

0.02536 

0.02262 

0.01961 

0.01630 






- 0.00883 

- 0.01088 

- 0.01209 

- 0.01262 

- 0.01266 

- 0.01232 

- 0.01169 

- 0.01085 

- 0.00992 

- 0.00896 

- 0.00807 

- 0.00728 

- 0.00662 

- 0.00607 

- 0.00561 

- 0.00528 

- 0.00498 

- 0.00466 

- 0.00425 

- 0.00370 

- 0.00298 

- 0.00204 

- 0.00085 

0 . 00060 

0.00233 

0.00436 

0.00669 

0.00933 

0.01227 

0.01550 






6.299 

2.939 

24.431 

0.0 

1 .0 



- 0.02043 

- 0.01732 

- 0.01391 

- 0.01025 

- 0.00643 

- 0.00254 

0.00136 

0.00526 

0.00914 

0.01299 

0.01679 

0.02052 

0.02412 

0.02756 

0.03079 

0.03376 

0.03643 

0.03875 

0.04065 

0.04208 

0.04301 

0.04340 

0.04321 

0.04247 

0.04122 

0.03950 

0.03737 

0.03484 

0.03196 

0.02874 






- 0.02043 

- 0.02293 

- 0.02453 

- 0.02531 

- 0.02538 

- 0.02486 

- 0.02389 

- 0.02251 

- 0.02081 

- 0.01885 

- 0.01670 

- 0.01443 

- 0.01212 

- 0.00983 

- 0.00760 

- 0.00542 

- 0.00332 

- 0.00124 

0.00082 

0.00290 

0.00501 

0 . 007 1 1 

0.00923 

0.01 139 

0.01366 

0.01608 

0.01870 

0.02155 

0.02463 

0.02794 






9.417 

4.409 

24.729 

0.0 

1 .0 



- 0.02314 

- 0.01979 

- 0.01627 

- 0.01259 

- 0.00879 

- 0.00491 

- 0.00098 

0.00295 

0.00690 

0.01087 

0.01488 

0.01895 

0.02307 

0.02723 

0.03139 

0.03549 

0.03945 

0.04319 

0.04655 

0.04948 

0.05189 

0.05375 

0.05503 

0.05571 

0.05582 

0.05540 

0.05446 

0.05300 

0.05105 

0.04864 






- 0.02314 

- 0.02593 

- 0.02782 

- 0.02884 

- 0.02907 

- 0.02861 

- 0.02757 

- 0.02605 

- 0.02412 

- 0.02185 

- 0.01930 

- 0.01650 

- 0.01350 

- 0.01035 

- 0.00708 

- 0.00373 

- 0.00031 

0.00319 

0.00673 

0.01030 

0.01389 

0.01747 

0.02105 

0.02464 

0.02826 

0.03198 

0.03580 

0.03970 

0.04372 

0.04784 






12.396 

5.879 

25.047 

0.0 

1 .0 



- 0.02485 

- 0.02140 

- 0.01783 

- 0.01417 

- 0.01039 

- 0.00654 

- 0.00263 

0.00135 

0.00537 

0.00944 

0.01359 

0.01784 

0.02220 

0.02669 

0.03127 

0.03590 

0.04048 

0.04491 

0.04907 

0.05288 

0.05626 

0.05915 

0.06154 

0.06344 

0.06486 

0.06576 

0.06618 

0.06608 

0.06551 

0.06447 






- 0.02485 

- 0.02779 

- 0.02982 

- 0.03097 

- 0.03131 

- 0.03091 

- 0.02985 

- 0.02824 

- 0.02616 

- 0.02370 

- 0.02091 

- 0.01784 

- 0.01452 

- 0.01098 

- 0.00724 

- 0.00333 

0.00072 

0.00491 

0.00925 

0.01371 

0.01826 

0.02287 

0.02756 

0.03237 

0.03729 

0.04234 

0.04751 

0.05279 

0.05817 

0.06367 






1 5 . 067 

7.349 

25.460 

0.0 

1 .0 



- 0.02641 

- 0.02293 

- 0.01929 

- 0.01551 

- 0.01158 

- 0.00751 

- 0.00333 

0.00108 

0.00560 

0.01014 

0.01458 

0.01891 

0.02332 

0.02781 

0.03239 

0.03702 

0.04164 

0.04616 

0.05049 

0.05459 

0.05841 

0.06196 

0.06518 

0.06804 

0.07050 

0.07253 

0.07412 

0.07525 

0.07591 

0.07612 






- 0.02641 

- 0.02935 

- 0.03133 

- 0.03240 

- 0.03259 

- 0.03197 

- 0.03065 

- 0.02859 

- 0.02600 

- 0.02306 

- 0.01996 

- 0.01680 

- 0.01343 

- 0.00988 
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- 0.00613 

- 0.00222 

0.00188 

0 

0.02568 

0.03120 

0.03696 

0 

0.06858 

0.07532 



17.461 

8.818 

26. 101 


- 0.02715 

- 0.02362 

- 0.01988 

-o 

0.00146 

0.00601 

0.01055 

0 

0.03204 

0.03639 

0.04075 

0 

0.06145 

0.06504 

0.06836 

0 

0.07979 

0 . 08090 



- 0.02715 

- 0.03004 

- 0.03192 

-0 

- 0.02821 

- 0.02559 

- 0.02264 

-0 

- 0.00648 

- 0.00285 

0.00099 

0 

0.02516 

0.03106 

0.03728 

0 

0.07245 

0.08010 



19.746 

10.288 

26.956 


- 0.02700 

- 0.02352 

- 0.01985 

-0 

0.00044 

0.00479 

0 . 009 1 7 

0 

0.03076 

0.03503 

0.03936 

0 

0.05973 

0.06331 

0.06665 

0 

0.07905 

0.08055 



- 0.02700 

- 0.02994 

- 0.03189 

-o 

- 0.02923 

- 0.02681 

- 0.02402 

-0 

- 0.00776 

- 0.00420 

- 0.00040 

0 

0.02345 

0.02933 

0.03557 

0 

0.07171 

0.07975 



22.012 

1 1 .758 

27.895 


- 0.03304 

- 0.02957 

- 0.02596 

-0 

- 0.00599 

- 0.00164 

0.00283 

0 

0.02572 

0.03008 

0.03435 

0 

0.05397 

0.05735 

0.06047 

0 

0.07205 

0.07348 



- 0.03304 

- 0.03599 

- 0.03800 

-0 

- 0.03566 

- 0.03324 

- 0.03037 

-0 

- 0.01280 

- 0.00916 

- 0.00541 

-0 

0.01769 

0.02337 

0.02939 

0 

0.06472 

0.07268 



24.277 

13.227 

28.850 


- 0.04749 

- 0.04402 

- 0.04047 

-0 

- 0.02092 

- 0.01662 

- 0.01219 

-0 

0.01 149 

0.01625 

0.02087 

0 

0 . 04007 

0.04320 

0.04603 

0 

0.05593 

0.05702 



- 0.04749 

- 0.05044 

- 0.05251 

-0 

- 0.05058 

- 0.04822 

- 0.04539 

-0 

- 0.02703 

- 0.02298 

- 0.01889 

-0 

0.00379 

0.00922 

0.01496 

0 

0.04860 

0.05622 



26.559 

14.697 

29 . 807 


- 0.07392 

- 0.07045 

- 0.06699 

-o 

- 0.04845 

- 0.04439 

- 0.04017 

-0 

- 0.01651 

- 0.01135 

- 0.00605 

-0 

0.01502 

0.01782 

0.02033 

0 

0.02809 

0.02867 



- 0.07392 

- 0.07687 

- 0.07903 

-0 

- 0.07814 

- 0.07599 

- 0.07337 

-0 

- 0.05503 

- 0.05058 

- 0.04581 

-0 

- 0.02128 

- 0.01616 

- 0.01076 

-0 

0.02076 

0.02787 




00616 

0.01067 

0.01541 

0.02040 

04294 

0.04911 

0.05546 

0.06195 

0.0 

1 .0 



01593 

- 0.01178 

- 0.00749 

- 0.00306 

01501 

0.01933 

0.02353 

0.02775 

04510 

0.04941 

0.05359 

0.05761 

07137 

0.07404 

0.07634 

0.07827 

03282 

- 0.03279 

- 0.03195 

- 0.03037 

01953 

- 0.01638 

- 0.01322 

- 0.00993 

00511 

0.00958 

0.01441 

0.01961 

04380 

0.05062 

0.05767 

0.06497 

0.0 

1 .0 



01605 

- 0.01210 

- 0.00802 

- 0.00384 

01359 

0.01800 

0.02233 

0.02655 

04366 

0.04788 

0.05199 

0.05596 

06974 

0.07253 

0.07504 

0.07720 

03294 

- 0.03311 

- 0.03248 

- 0.03116 

02095 

- 0.01771 

- 0.01442 

- 0.01 113 

00366 

0.00806 

0.01282 

0.01795 

04217 

0.0491 1 

0.05637 

0.06391 

0.0 

1 .0 



02222 

- 0.01835 

- 0.01436 

- 0.01025 

00738 

0.01 199 

0.01663 

0.02123 

03856 

0.04267 

0.04661 

0.05039 

06335 

0.06595 

0.06828 

0.07031 

03911 

- 0.03936 

- 0.03882 

- 0.03757 

02717 

- 0.02372 

- 0.0201 1 

- 0.01645 

00143 

0.00284 

0.00743 

0.01238 

03579 

0.04253 

0.04961 

0.05701 

0.0 

1 .0 



03680 

- 0.03302 

- 0.0291 1 

- 0.02508 

00763 

- 0.00294 

0.00185 

0.00667 

02526 

0.02930 

0.03313 

0.03672 

04860 

0.05089 

0.05286 

0.05455 

05368 

- 0.05403 

- 0.05357 

- 0.05239 

04218 

- 0.03865 

- 0.03490 

- 0.03102 

01473 

- 0.01053 

- 0.00605 

- 0.00128 

02103 

0.02746 

0.03420 

0.04125 

0.0 

1.0 



06346 

- 0.05982 

- 0.05617 

- 0.05236 

03584 

- 0.03127 

- 0.02650 

- 0.02161 

00091 

0.00404 

0.00839 

0.01188 

02252 

0.02443 

0.02596 

0.02716 

08034 

- 0.08081 

- 0.08063 

- 0.07969 

07036 

- 0.06698 

- 0.06325 

- 0.05928 

04093 

- 0.03579 

- 0.03077 

- 0.02614 

00502 

0.00103 

0.00730 

0.01386 



OUTPUT DATA FORMAT 

Output data format for the TSCLP Code is very similar to that for the basic 
NASA/Grumman Transonic Wing-Body Codes (Refs. 1,2). Geometry verification is ob- 
tained only by specifying zero iterations. Printed and plotted output includes 
information sufficient to verify input configuration geometry, generated grid struc- 
tures, and computational representation of configuration components. Analysis runs 
are obtained when some non-zero number of iterations is specified. Printed output 
includes solution iteration convergence history followed by resultant surface veloc- 
ities and pressures, load distributions, and force and moment coefficients for each 
component. Plotted output of surface pressures, load distributions, and force and 
moment coefficients is produced only when some non-zero number of "fine grid" 
iterations is specified. Printed and plotted output will then also include a force 
and moment summary for the entire configuration, on a component-by-component basis. 

Note that integrated forces and moments are output both with and without store 
body viscous crossflow estimates. Viscous skin friction estimates are not included 
in the integrated forces and moments but are always output as separate, distinct 
quantities. 
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SUBROUTINE CALL SEQUENCE 


MAIN PROGRAM SUBROUTINE CALL SEQUENCE : STATS 


rsTarn 


WNGDAT 

— BODDAT r -rWlKDE^ 

U~owrmr^ 

- CRUSET 1 — BODSET — 


U WNGSET 


— STRDAT 

— WSISET 

— CEESET 


— BUMP 


M QWIKLQ 
SLOPE 


U SPLINE 

L [qwTklo"]- > 


- SLOPE - 
t- SNGLRT 

- CEEDIF 
— ARC2D 

- INT2D 
— INT3D 
— FINTIP 


U PYLSET 


— (contrlI— » 

— PLTSET — 


> — FINSET 

E lgrng 
PYLBCS 
CBCSPL 

— [qwiklo L -4 


SLOPE 

— BODLIM 

- SAREA 

— j QW I KL Q > 


SPLINE 


SPLINE 

INT3D 


E 


BUMP 

LGRNG 

FINBCS 


- SPLINE 
f QWiKLO~i -> 


CBCSPL 
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SUBROUTINE CALL SEQUENCE : CONTRL 


- 4 contrlT 


OUTPUT 


CRUOUT 

— WSIOUT 

— WFGOUT — 

— SMOTH 

— 1 STRIPK 1 — ^ 

— BODOUT 


— GLOBAL 

— SUPSET 

— WINGF - 


[- I^WFGFIL T -» 
BODVAL 
BLAR 


r 


U WSIREL 


— CE1REL 


t 


PYLCPS 

CE1CPS 

CE2CPS 

FMSOUT 

CRUCUT 

TRID 

MCHANG 

TRID 


SPLINE 
SMOTH 
VINTER 
STRIPK j— » 

WSICUT 

TRID 

CE1PYL 


U CE2REL 


ARC1 

ARC 2 

ARC3 

TRIDBC 

CE2PYL 

ARC2 

ARC3 

TRIDBC 


— CRUINT 

— WSIINT 

— CE1INT 

- 4 CRLFIL~1 — » 
WSIFIL ) -> 
- I CEIFIL U 

— CE2FIL 

— FINFIL 

— PYLFIL 


\ WFGFIL } -» 


SPLINE 

BODFM 

QWIKLO l— > 

CE1REL 

TRIDBC 

CE2REL 

TRIDBC 

TRID 



— TRID 

j- ARC1 
I— ARC 2 


ARC 2 



SUBROUTINE CALL SEQUENCE : WFGFIL 


SUBROUTINE CALL SEQUENCE : CRUFIL 



WFGP 

WFGV 

WFGVX 

WFGVZ 

WFGZ2 

WFGZ3 

WFGZ4 

WFGZ5 

WFGZ6 

WFGZ7 

CRUZ2 

CRUZ3 

CRUZ4 

CRUZ 5 

CRUZ6 

CRUZ7 

WSIZ2 

WSIZ3 

WSIZ4 

WSIZ5 

WSIZ6 

WSIZ7 

TRIDBC 



TRIDBC 

TRIDBC 

TRIDBC 

TRIDBC 



CRUZ 2 
CRUZ 3 
CRUZ 4 
CRUZ 5 
CRUZ 6 
CRUZ 7 
WSIZ2 
WSIZ3 
WSIZ4 
WSIZ5 
WSIZ6 
WSIZ7 


SUBROUTINE CALL SEQUENCE : WSIFIL SUBROUTINE CALL SEQUENCE : CE1FIL 



CRUZ 2 
CRUZ3 
CRUZ 4 
CRUZ 5 
CRUZ6 
CRUZ 7 
WSIZ2 
WSIZ3 
WSIZ4 
WSIZ5 
WSIZ6 
WSIZ7 
CE1Z4 
CE1Z5 
CE1Z6 
CE1Z7 



WSIZ2 

WSIZ3 

WSIZ4 

WSIZ5 

WSIZ6 

WSIZ7 

CE1Z4 

CE1Z5 

CE1Z6 

CE1Z7 
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SUBROUTINE CALL SEQUENCE : QWIKDE 


-1 OWIKDE 


CSMDEF DSETUP 

I— DLOKUP 


— CSMCHK 


— BLMDEF 


DSETUP 

DLOKUP 

CURVES 

KRVDEF 


— BLMCHK 
I— GEMOUT 


SUBROUTINE CALL SEQUENCE : QWIKLO 


—{QWIKLO 


CSGEOM 


L_ SMOTH 


BLMSET 
BLGEOM 
CSMSET 

U THELIM 


H 


CURVES 

BLGEOM 

CSMCOE 

AORDER 

CSMFLT 


I— CSMINT 


CSCALC 


BLGEOM 

CSMCOE 

SINCOS 

LINLIN 

LINELL 

ELLELL 


L- CSCALC 


E SINCOS 
VDOTV 
MDOTV 


E SINCOS 
VDOTV 
MDOTV 


— CROSS 
y- ELLCAL 
L SINCOS 
j- ELLCAL 
— CROSS 
— ESTNXT 
— SETNXT 


SINCOS 

SINCOS 
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SUBROUTINE CALL SEQUENCE : STRIPK 


STRIPK I- 


BLLAM 


I — VBRAD 


E SLOPY ■ 
SPLN1 
ORDIN 
y- VNUSUB 
— SLOPBL 
- TANCAL 
— SPRINT 
— FINT 
— SOLVEB 
— ORDIN 
- REDUCX 
— RLORD 
— GORD 
- SIMPSN 


SPLN1 


GRAD 

GORD 



PLOTTING PROGRAM SUBROUTINE CALL SEQUENCE : STATSP 


t STATSP l 


GEVER 

STRVER 

WFGPLT 

BODPLT 

PYLPLT 

CE2PLT 

FMSPLT 


HIGH-LEVEL PLOTTING SUBROUTINES CALLED BY THE ABOVE 
CRAY XMP versions contain DISSPLA calls 
LaRC CDC versions contain Langley CALCOMP calls 


PGSTRT 

PGSTOP 

PSTRT 

PSTOP 

SPSTRT 

SPSTOP 

GSETUP 

XSETUP 

YSETUP 

SETHIT 

WMESS 

WREAL 

WINTG 

DRAWP 

DRAWG 



SUBROUTINE DESCRIPTION 


AORDER 

ARC1 

ARC2 

ARC2D 

ARC3 


Orders a set of numbers by permutation index 
Inverts tridiagonal matrix for C-grids 

Inverts tridiagonal matrix for C-grids (modified storage scheme) 

Computes arc length along a 2-D curve 

Inverts periodic tridiagonal matrix for C-grids 


BLAR 


BLGEOM 

BLUM 

BLMCHK 

BLMDEF 

BLMSET 

BLOCK DATA 

BODDAT 

BODFM 

BODLIM 

BODOUT 

BODPLT 

BODSET 

BODVAL 

BUMP 


Main control routine for laminar and turbulent modified chordwise 
boundary layer calculation. Computes boundary layer displacement 
thickness slope for vi scous/invi scid interaction mode of operation 
Assigns body line model values and derivatives to control point coor- 
dinates 

Computes Thwaites laminar boundary ‘layer with Rott and Crabtree 
compre s s i bi 1 i ty modification 

Correlates and checks the input data deck and the indices for the 

generated body line math models 

Defines body line models from the input data 

Controls the determination values and first and second derivatives 

for all body line models at a given x-station 

Initializes in-core memory to zero 

Reads axi symmetric or QUICK fuselage input data 

Computes integrated fuselage force and moment coefficients 

Computes fuselage global crude grid limiters 

Prints out solution for fuselage 

Plots solution for fuselage 

Controls set up for fuselage in global crude grid 
Computes fuselage boundary point potential values 
Bump function for grid modifications 


CBCSPL 

CE1CPS 

CE1FIL 

CE1INT 

CE1PYL 

CE1REL 


Computes cubic splines through sets of points 
Prints out coarse C-grid solution for store 
Interpolation routine for coarse C-grid boundaries 
Initializes interpolation parameters for use in CE1FIL 
Relaxation routine to treat pylon surface in coarse C-grid 
Relaxation routine for coarse C-grid 
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CE1Z4 

CE1Z5 

CE1Z6 

CE1Z7 

CE2CPS 

CE2FIL 

CE2PLT 

CE2PYL 

CE2REL 

CEEDIF 

CEESET 

CONTRL 

CROSS 

CRUCUT 

CRUFIL 

CRUINT 

CRUOUT 

CRUSET 

CRUZ2 

CRUZ3 

CRUZ4 

CRUZS 

CRUZ6 

CRUZ7 

CSCALC 


Zone 4 set up for coarse C-grid interpolations (above wing and out- 
board of pylon) 

Zone 5 set up for coarse C-grid interpolations (above wing and 
inboard of pylon) 

Zone 6 set up for coarse C-grid interpolations (below wing and out- 
board of pylon) 

Zone 7 set up for coarse C-grid interpolations (below wing and 
inboard of pylon) 

Prints out fine C-grid solution for store 

Interpolation routine for fine C-grid boundaries 

Plots fine C-grid solution for store 

Relaxation routine to treat pylon surface in fine C-grid 

Relaxation routine for fine C-grid 

Computes C-grid metrics 

Sets up cylindrical C-grids 

Main control routine for relaxation solution, interpolations, bound- 
ary layer analysi s, and printed output 
Solves for the intersection of two lines in a plane 
Relaxation routine to treat global coarse grid inner boundary 
Interpolation routine for global coarse grid inner boundary 
Initializes interpolation parameters for use in CRUFIL 
Prints out global coarse grid solution for wing 
Set up global coarse grid 

Zone 2 set up for global coarse grid interpolations (above wing) 

Zone 3 set up for global coarse grid interpolations (below wing) 

Zone 4 set up for global coarse grid interpolations (above wing and 
outboard of pylon) 

Zone 5 set up for global coarse grid interpolations (above wing and 
inboard of pylon) 

Zone 6 set up for global coarse grid interpolations (below wing and 
outboard of pylon) 

Zone 7 set up for global coarse grid interpolations (below wing and 
inboard of pylon) 

Computes radial position and derivatives for specified cross-section 
model 
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CSGEOM 

CSMCHK 

CSMCOE 

CSMDEF 

CSMFLT 

CSMINT 

CSMSET 

CURVES 


Is the main subroutine in the look-up portion of the QUICK system. 
It calls appropriate subroutines to evaluate body line values and 
construct cross-section geometry at a given x-station. It is used 
for all geometry model interrogation 

Correlates and checks the input data deck and the indices for the 
cross-sectional math model 

Composes the equations which are to define the cross-section geometry 
at a given x-station 

Logically defines the cross-section models from the input data 
Creates control point definitions to permit the insertion of a smooth 
fillet between cross-sectional arcs 

Locates user specified intersections between cross-sectional arcs and 
adjusts their use (theta limits) 

Sets up the control point coordinate arrays used to define the 
cross-section geometry at a specified x-station 

Calculates values and first and second derivatives for individual 
curve fits 


DLOKUP 


DRAWG 

DRAWP 

DSETUP 


ELLCAL 

ELLELL 

ESTNXT 

FINBCS 

FINFIL 

FINSET 

FINT 


Is a simple dictionary look-up routine. It assigns an index to match 
an input name to a codeword list, but is not capable of adding new 
i terns to that 1 i st 

Plots a curve according to graph axes 

Plots a curve according to plot area dimensions 

Is an adapting dictionary look-up routine. New items are added to a 
codeword list, an index is returned for the codeword, and an indica- 
tor (INEW) is set equal to 1 when a new item is encountered 

Set up for el 1 ipse 

Calculates intersection of two ellipses 

Estimates non-linear root by modified inverse quadratic 

Sets up store fin surface boundary conditions in C-grids 
Interpolation routine for store fin surface potentials in coarse 
C-grid 

Sets up store fin geometry in C-grids 
Simultaneous triple interpolation 


102 



FINTIP 

FMSOUT 

FMSPLT 

GEMOUT 

GEVER 

GLOBAL 

GORD 

GRAD 

GSETUP 

INT2D 

INT3D 

KRVDEF 


LGRNG 

LINELL 

LINLIN 

MCHANG 

MDOTV 

ORDIN 

OUTPUT 

PGSTOP 

PGSTRT 

PLTSET 

PSTOP 

PSTRT 

PYLBCS 

PYLCPS 

PYLFIL 


Modifies C-grids to follow store fin tip vortex streamlines 
Prints out summary of configuration force and moment coefficients 
Plots summary of configuration force and moment coefficients 

Ensures that all body lines required by a cross-sectional model are 
defined for the range of that model 

Plots geometry verification for wing, fuselage, and pylon 
Relaxation routine for global coarse grid 
Bradshaw' s G function 

Slope of a function at its tabulated points 
Sets up graph area size and axis scales 

Interpolates along a 2-D curve 
Interpolates along a 3-D curve 

Calculates coefficients for the various curve fits associated with 
body line math models 

Interpolates using Lagrange polynomials 

Solves for the intersection of a line and an ellipse 

Solves for the intersection of two lines 

Computes local Mach angle from coefficients of flow equation 
Performs matrix multiplication of a vector 

Linear interpolation 
Controls print out of solution 

Terminates plotting system 

Initiates plotting system 

Set up routine for plotting program 

Terminates a plot 

Initiates a plot 

Sets up pylon surface boundary conditions 
Prints out solution for pylon 

Interpolation routine for pylon surface potentials in coarse C-grid 
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PYLPLT 

PYLSET 

QWIKDE 

QWIKLO 

REDUCX 

RLORD 

SAREA 

SETHIT 

SETNXT 

SIMPSN 

SINCOS 

SLOPBL 

SLOPE 

SLOPY 

SMOTH 

SNGLRT 

SOLVEB 

SPLINE 

SPLN1 

SPRINT 

SPSTOP 

SPSTRT 

STATS 

STATSP 

STRDAT 

STRIPK 

STRVER 

SUPSET 


Plots solution for pylon 
Sets up pylon geometry 

Main control routine for Quick geometry definition and check out 
Main control routine for interrogation of Quick geometry math model 

Performs interpolation to new grid 
Bradshaw's L function 

Computes body surface area given an array of cross-sections 
Sets height for text appearing on plots 
Reorders points for non-linear root finder 
Simpson's integration rule 

Adjusts input interrogation angles for top and bottom dead center 

Slope of a tabulated function at an arbitrary point 

Computes boundary conditions for wing surface and axi symmetric bodies 

Computes wing surface slopes 

Function for smoothing an array of values 

Used for internal calculation of singularity location in C-grid 
conformal mapping 

Solution of two simultaneous linear algebraic equations 
Computes a cubic spl ine through a set of points 

Computes continuous derivatives interpolation by means of a cubic fit 

Prints output of profile results 

Terminates a subplot 

Initiates a subplot 

Main program 

Main plotting program 

Reads store and pylon input data 

Starting condition setup and flow control for laminar/turbulent 
boundary layer prediction 

Plots geometry and grid set up verification for stores 

Applies supersonic inflow, outflow, and radiation-type boundary 

conditions 
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TANCAL 

THELIM 

TRID 

TRIDBC 

VBRAD 


VDOTV 

VINTER 

VNUSUB 

WFGFIL 

WFGOUT 

WFGP 

WFGPLT 

WFGV 

WFGVX 

WFGVZ 

WFGZ2 

WFGZ3 

WFGZ4 

WFGZ5 

WFGZ6 

WFGZ7 

WINGF 

WINTG 


Computes characteristic angles for use in boundary layer solution 
Creates and controls use of theta arrays to establish continuity in 
the cross-section model 
Solves tridiagonal matrix 

Solves tridiagonal matrix with special end conditions 

Computes Bradshaw's compressible 2-D turbulent boundary layer 
simulating 3-D boundary layer on infinite yawed wing by Nash-Tseng 
modified chord technique 
Computes a vector dot product 

Performs cubic fit for separated boundary layer in wing section cove 
regions 

Computes the Nash effective viscosity 

Interpolation routine for wing fine grid boundaries 
Prints out wing fine grid solution 
Dirichlet-type boundary conditions for wing fine grid 
Plots wing fine grid solution 

Neumann-type boundary conditions for wing fine grid cut-out due to 
store 

Neumann-type boundary conditions for wing fine grid upstream and 
downstream boundaries 

Neumann-type boundary conditions for wing fine grid top and bottom 
boundaries 

Zone 2 set up for wing fine grid interpolations (above wing) 

Zone 3 set up for wing fine grid interpolations (below wing) 

Zone 4 set up for wing fine grid interpolations (above wing and 

outboard of pylon) 

Zone 5 set up for wing fine grid interpolations (above wing and 

inboard of pylon) 

Zone 6 set up for wing fine grid interpolations (below wing and 

outboard of pylon) 

Zone 7 set up for wing fine grid interpolations (below wing and 

i nboard of pyl on) 

Relaxation routine for wing fine grid 
Writes an integer number on a plot. 
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WMESS Writes a character string on a plot 

WNGDAT Reads input data for wing 

WNGSET Sets up wing geometry 

WREAL Writes a real number on a plot 

WSICUT Relaxation routine to treat WSI grid inner boundary 

WS I FI L Interpolation routine for WSI grid inner and outer boundaries 

WSI INT Initializes interpolation parameters for use in WSI FI L 

WSIOUT Prints out WSI grid solution for wing 

WSIREL Relaxation routine for WSI grid 

WSISET Sets up WSI grid 

WSIZ2 Zone 2 set up for WSI grid interpolations (above wing) 

WSIZ3 Zone 3 set up for WSI grid interpolations (below wing) 

WSIZ4 Zone 4 set up for WSI grid interpolations (above wing and outboard of 

pylon) 

WSIZ5 Zone 5 set up for WSI grid interpolations (above wing and inboard of 

pylon) 

WSIZ6 Zone 6 set up for WSI grid interpolations (below wing and outboard of 

pylon) 

WSIZ7 Zone 7 set up for WSI grid interpolations (below wing and inboard of 

pylon) 

XSETUP Draws and labels x-axis 


YSETUP Draws and labels y-axis 
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KEY VARIABLE DESCRIPTION 


AAXIS1 

AAXIS2 

AK 

ALPAS 

ALPHA 

AM2 

AMACH 

ANGF 

AOA 

BAREA 

BAX IS 

BAX I SO 

BCF 

BCL 

BCLB 

BCLF 

BCU 

BCUB 

BCUF 

BDD 

BETAP 

BETAS 

BLD 

BNOSE 

BODCD 

BODCL 

BODCM 

BPAREA 

BS 

BTAIL 


Nominal extent of C-grids upstream of store nose, as fraction of store 
body length 

Nominal extent of C-grids downstream of store tail, as fraction of 
store body length 
The value 1-M^ 

Store pitch angle relative to aircraft 

Angle-of-attack (radians) 

2 

The value M 

Freestream Mach number 
Store fin angular locations 
Angle-of-attack (degrees) 

Fuselage wetted area 
Nominal radius of C-grids 

Nominal radius inherent in C-grid conformal mapping 

Integrated fuselage skin friction coefficient 

Wing global coarse grid lower boundary slopes 

Wing WSI grid lower boundary slopes 

Wing fine grid lower boundary slopes 

Wing global coarse grid upper boundary slopes 

Wing WSI grid upper boundary slopes 

Wing fine grid upper boundary slopes 

Fuselage longitudinal distribution of axial force 

Pylon yaw angle relative to aircraft 

Store yaw angle relative to aircraft 

Fuselage longitudinal distribution of normal force 

X coordinate of fuselage nose 

Integrated fuselage drag coefficient 

Integrated fuselage lift coefficient 

Integrated fuselage pitching moment coefficient 

Fuselage projected planform area 

Fuselage plot scaling parameter 

X coordinate of fuselage tail 
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CA 

Global coarse grid metric £ x 


CAB 

WSI grid metric £ x 


CAV 

Wing average chord, c a ^ Q 


CB 

Global coarse grid metric £ xx 


CBB 

WSI grid metric £ 

XX 


CC 

Global coarse grid metric ti 

y 


CCB 

WSI grid metric n 

y 


CD 

Global coarse grid metric ti 

yy 


CDB 

WSI grid metric ri 

yy 


CDINT 

Integrated wing section drag coefficient, c^ 


CE 

Global coarse grid metric C z 


CEB 

WSI grid metric c z 


CF 

Global coarse grid metric c 


CFB 

WSI grid metric c zz 


CFINT 

Integrated wing section skin friction coefficient, c^ 


CIR 

Wing circulation in global coarse and wing fine grids 


CIRB 

Wing circulation in WSI grid 


CIRFC 

Fin circulation in coarse C-grid 


CIRFS 

Fin circulation in fine C-grid 


CLINT 

Integrated wing section lift coefficient, c^ 


CMINT 

Integrated wing section pitching moment coefficient, c m about REFX 


CMLOC 

Integrated wing section pitching moment coefficient 


CMXFIN 

X component of store fin moment coefficient (body axes) 


CMXFNV 

X component of store fin moment coefficient (body axes) including 
friction and viscous crossflow 

ski n 

CMXSBD 

X component of store body moment coefficient (body axes) 


CMXSBV 

X component of store body moment coefficient (body axes) including skin 
friction and viscous crossflow 

CMXSTR 

X component of store moment coefficient (body axes) 


CMXSTV 

X component of store moment coefficient (body axes) including 
friction and viscous crossflow 

skin 

CMYFIN 

Y component of store fin moment coefficient (body axes) 


CMYFNV 

Y component of store fin moment coefficient (body axes) including 
friction and viscous crossflow 

ski n 

CMYSBD 

Y component of store body moment coefficient (body axes) 
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CMYSBV 

CMYSTR 

CMYSTV 

CMZFIN 

CMZFNV 

CMZSBD 

CMZSBV 

CMZSTR 

CMZSTV 

COST 

CPL 

CPU 

CSCUT 

CXFIN 

CXFNV 

CXSBD 

CXSBV 

CXSTR 

CXSTV 

CYFIN 

CYFNV 

CYSBD 

CYSBV 

CYSTR 

CYSTV 


Y component of store body moment coefficient (body axes) including skin 
friction and viscous crossflow 

Y component of store moment coefficient (body axes) 

Y component of store moment coefficient (body axes) including skin 
friction and viscous crossflow 

Z component of store fin moment coefficient (body axes) 

Z component of store fin moment coefficient (body axes) including skin 
friction and viscous crossflow 

Z component of store body moment coefficient (body axes) 

Z component of store body moment coefficient (body axes) including skin 
friction and viscous crossflow 
Z component of store moment coefficient (body axes) 

Z component of store moment coefficient (body axes) including skin 

friction and viscous crossflow 

Trigonometric cosine of C-grid 0 coordinates 

Lower surface pressure coefficient 

Upper surface pressure coefficient 

Fuselage X station for cross-section cut 

X component of store fin force coefficient (body axes) 

X component of store fin force coefficient (body axes) including skin 
friction and viscous crossflow 

X component of store body force coefficient (body axes) 

X component of store body force coefficient (body axes) including skin 
friction and viscous crossflow 
X component of store force coefficient (body axes) 

X component of store force coefficient (body axes) including skin fric- 
tion and viscous crossflow 

Y component of store fin force coefficient (body axes) 

Y component of store fin force coefficient (body axes) including skin 
friction and viscous crossflow 

Y component of store body force coefficient (body axes) 

Y component of store body force coefficient (body axes) including skin 
friction and viscous crossflow 

Y component of store force coefficient (body axes) 

Y component of store force coefficient (body axes) including skin fric- 
tion and viscous crossflow 
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CZFIN 

CZFNV 

CZSBD 

CZSBV 

CZSTR 

CZSTV 


DCYDX 

DCYDXV 

DCZDX 

DCZDXV 

DELFL 

DELFP 

DELFU 

DELGL 

DELGP 

DELGU 

DELML 

DELMP 

DELMU 

DELTAF 

DETA 

DETAB 

DETAS 

DETASS 

DIM 

DRDXC 

DTS 


Z component of store fin force coefficient (body axes) 

Z component of store fin force coefficient (body axes) including skin 
friction and viscous crossflow 

Z component of store body force coefficient (body axes) 

Z component of store body force coefficient (body axes) including skin 
friction and viscous crossflow 
Z component of store force coefficient (body axes) 

Z component of store force coefficient (body axes) including skin fric- 
tion and viscous crossflow 

Store body longitudinal distribution of side force 

Store body longitudinal distribution of side force, including viscous 
crossflow 

Store body longitudinal distribution of normal force 

Store body longitudinal distribution of normal force, including viscous 
crossflow 

Wing lower surface boundary layer slopes in wing fine grid 

Wing surface boundary layer slopes in wing fine grid at inboard 

wing/pylon junction 

Wing upper surface boundary layer slopes in wing fine grid 

Wing lower surface boundary layer slopes in global coarse grid 

Wing surface boundary layer slopes in global coarse grid at inboard 
wing/pylon junction 

Wing upper surface boundary layer slopes in global coarse grid 

Wing lower surface boundary layer slopes in WSI grid 

Wing surface boundary layer slopes in WSI grid at inboard wing/pylon 
junction 

Wing upper surface boundary layer slopes in WSI grid 

Individual, al 1-moveabl e , store fin deflection angle 

Global coarse grid n mesh cell size 

WSI grid n mesh cell size 

Coarse C-grid n mesh cell size 

Fine C-grid n mesh cell size 

Reference length for non-dimensional izing maximum change in potential 
Axisymmetric fuselage slope distribution in global coarse grid 
C-grid 0 mesh cell size 
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DXI 

Global coarse grid £ mesh cell size 


DXIB 

WSI grid £ mesh cell size 


DXSIS 

Coarse C-grid £ mesh cell size 


DXSISS 

Fine C-grid £ mesh cell size 


DXW 

Wing fine grid X mesh cell size 


DZETA 

Global coarse grid c mesh cell size 


DZETAB 

WSI grid £ mesh cell size 


DZW 

Wing fine grid Z mesh cell size 


ETA 

Global coarse grid n coordinates 


ETAB 

WSI grid x\ coordinates 


ETRRS 

Coarse C-grid metric ri rr 


ETRRSS 

Fine C-grid metric n 


ETRS 

Coarse C-grid metric n r 


ETRSS 

Fine C-grid metric n r 


ETXRS 

Coarse C-grid metric ti 

a xr 


ETXRSS 

Fine C-grid metric n xr 


ETXS 

Coarse C-grid metric n x 


ETXSS 

Fine C-grid metric n x 


ETXXS 

Coarse C-grid metric ti xx 


ETXXSS 

Fine C-grid metric ti 

X X 


FINCMX 

X component of store fin moment coefficient (stability axes) 


FINCMY 

Y component of store fin moment coefficient (stability axes) 


FINCMZ 

Z component of store fin moment coefficient (stability axes) 


FINCX 

X component of store fin force coefficient (stability axes) 


FINCY 

Y component of store fin force coefficient (stability axes) 


FINCZ 

Z component of store fin force coefficient (stability axes) 


FNSCF 

Integrated store fin section skin friction coefficient 


FNSCMY 

Integrated store fin section pitching moment coefficient 


FNSCX 

Integrated store fin section axial force coefficient 


FNSCZ 

Integrated store fin section normal force coefficient 


FNVCMX 

X component of store fin moment coefficient (stability axes) 
skin friction and viscous crossflow 

i ncl udi ng 

FNVCMY 

Y component of store fin moment coefficient (stability axes) 
skin friction and viscous crossflow 

i ncl udi ng 
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FNVCMZ 

FNVCX 

FNVCY 

FNVCZ 

FNXC 

FNXS 

FNYC 

FNYS 

FRLEC 

FRLES 

FRTEC 

FRTES 

FS 

FXLEC 

FXLES 

FXTEC 

FXTES 

G 

H 

IBCL 

ICE10 

ICRUI 

IFLEC 

IFLES 

IFTEC 

IFTES 

IL 

I LB 

ILEF 


Z component of store fin moment coefficient (stability axes) including 
skin friction and viscous crossflow 

X component of store fin force coefficient (stability axes) including 
skin friction and viscous crossflow 

Y component of store fin force coefficient (stability axes) including 
skin friction and viscous crossflow 

Z component of store fin force coefficient (stability axes) including 
skin friction and viscous crossflow 
X component of store fin surface normal in coarse C-grid 

X component of store fin surface normal in fine C-grid 

R component of store fin surface normal in coarse C-grid 

R component of store fin surface normal in fine C-grid 

R coordinate of store fin section leading edge in coarse C-grid 

R coordinate of store fin section leading edge in fine C-grid 

R coordinate of store fin section trailing edge in coarse C-grid 

R coordinate of store fin section trailing edge in fine C-grid 

Store fin plot scaling parameter 

X coordinate of store fin section leading edge in coarse C-grid 

X coordinate of store fin section leading edge in fine C-grid 

X coordinate of store fin section trailing edge in coarse C-grid 

X coordinate of store fin section trailing edge in fine C-grid 

The val ue (2f+l)M^ 


The value (y-l)M‘ 


WSI grid I value closest to store nose 

I interpolation index to locate coarse C-grid point in WSI grid 
I interpolation index to locate global coarse grid point in WSI grid 
I value of store fin section leading edge point in coarse C-grid 

I value of store fin section leading edge point in fine C-grid 

I value of store fin section trailing edge point in coarse C-grid 

I value of store fin section trailing edge point in fine C-grid 

Global coarse grid wing leading edge I value 
WSI grid wing leading edge I value 
Wing fine grid wing leading edge I value 
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IMACH 

I MAX 

IMAXB 

IMAXS 

IMAXSS 

IMAXW 

INOSEC 

IPLEC 

IPLEF 

IP LEG 

IPLEM 

IPLES 

IPTEC 

IPTEF 

IPTEG 

IPTEM 

IPTES 

IT 

ITAILC 

ITB 

ITEF 

ITER 

IWFG 

IWSIF 

IWSII 

IWSIO 

JBCL 

JBLI 

JBLO 

JBRI 

JBRO 

JCEIO 


Code for subsonic (0) or supersonic (1) flow at a grid point 

Number of global coarse grid points in X direction 

Number of WSI grid points in X direction 

Number of coarse C-grid £ mesh cell points 

Number of fine C-grid £ mesh cell points 

Number of wing fine grid points in X direction 

Global coarse grid I value at fuselage nose 

I value of pylon section leading edge point in coarse C-grid 

I value of pylon section leading edge point in wing fine grid 

I value of pylon section leading edge point in global coarse grid 

I value of pylon section leading edge point in WSI grid 

I value of pylon section leading edge point in fine C-grid 

I value of pylon section trailing edge point in coarse C-grid 

I value of pylon section trailing edge point in wing fine grid 

I value of pylon section trailing edge point in global coarse grid 

I value of pylon section trailing edge point in WSI grid 

I value of pylon section trailing edge point in fine C-grid 

Global coarse grid wing trailing edge I value 

Global coarse grid I value at fuselage tail 

WSI grid wing trailing edge I value 

Wing fine grid wing trailing edge I value 

Iteration count 

I interpolation index to locate wing fine grid point in global coarse 
grid 

I interpolation index to locate WSI inner-boundary front-face grid 
point in coarse C-grid 

I interpolation index to locate WSI inner-boundary side-face grid point 
in coarse C-grid 

I interpolation index to locate WSI grid point in global coarse grid 

WSI grid J value at store centerline 

Inner J location of inboard WSI grid inner boundary 

Outer J location of inboard WSI grid inner boundary 

Inner J location of outboard WSI grid inner boundary 

Outer J location of outboard WSI grid inner boundary 

J interpolation index to locate coarse C-grid point in WSI grid 
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JCRUI 

JFINEI 

JFINEO 

JFINTC 

JFINTS 

JMAX 

JMAXB 

JMAXS 

JMAXSS 

JPYLTC 

JPYLTS 

JROOT 

JSDC 

JSTOR 

JSTORI 

JSTORO 

JTIP 

JWFG 

JWSI 

JWSIF 

JWSII 

JWSIO 

KBB 

KBBOT 

KBC 

KBCL 

KBTOP 

KBW 

KCEIO 

KCRUI 

KFINEB 

KFINET 


J interpolation index to locate global coarse grid point in WSI grid 

Inboard J location of wing fine grid cutout due to store 

Outboard J location of wing fine grid cutout due to store 

J value of store fin tip in coarse C-grid 

J value of store fin tip in fine C-grid 

Number of global coarse grid points in Y direction 

Number of WSI grid points in Y direction 

Number of coarse C-grid n mesh cell points 

Number of fine C-grid n mesh cell points 

J extent of pylon in coarse C-grid 

J extent of pylon in fine C-grid 

Global coarse grid J value at wing root 

Global coarse grid J value just outboard of wing/fuselage junction 
Global coarse grid J value of pylon and/or store centerline 
Inboard J location of global coarse grid inner boundary 
Outboard J location of global coarse grid inner boundary 
Global coarse grid J value at wing tip 

J interpolation index to locate wing fine grid point in global coarse 
grid 

J array to order WSI grid inner-boundary side-face J-K grid points 
J interpolation index to locate WSI inner-boundary front-face grid 
point in coarse C-grid 

J interpolation index to locate WSI inner-boundary side-face grid point 
in coarse C-grid 

J interpolation index to locate WSI grid point in global coarse grid 

WSI grid K value at wing plane 

Lower K location of WSI grid inner boundary 

Global coarse grid K value at wing plane 

WSI grid K value at store centerline 

Upper K location of WSI grid inner boundary 

Wing fine grid K value at wing plane 

K interpolation index to locate coarse C-grid point in WSI grid 
K interpolation index to locate global coarse grid point in WSI grid 
Lower K location of wing fine grid cutout due to store 
Upper K location of wing fine grid cutout due to store 
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KFS 

KLO 

KLOC 

KMAX 

KMAXB 

KMAXS 

KMAXW 

KODB 

KPS 

KPYLRF 

KPYLRG 

KPYLRM 

KPYLTF 

KPYLTG 

KPYLTM 

KREFL 

KSTORB 

KSTORT 

KUP 

KUPC 

KWFG 

KWSI 

KWSIF 

KWSI I 

KWSIO 

MAX IT 
MAXITF 
MAXITM 
MODV 

NANGFS 

NBODY 


K value of store fin in C-grid 
K indice of 6+270° coordinate 

Global coarse grid K limiters for lower portion of fuselage 

Number of global coarse grid points in Z direction 

Number of WSI grid points in Z direction 

Number of C-grid 0 mesh cell points 

Number of wing fine grid points in Z direction 

Axi symmetric or QUICK fuselage definition flag 

K value of pylon in C-grid 

Pylon lower K limit in wing fine grid 

Pylon lower K limit in global coarse grid 

Pylon lower K limit in WSI grid 

Pylon upper K limit in wing fine grid 

Pylon upper K limit in global coarse grid 

Pylon upper K limit in WSI grid 

K indice of 0+180° coordinate 

Lower K location of global coarse grid inner boundary 
Upper K location of global coarse grid inner boundary 
K indice of 0+90° coordinate 

Global coarse grid K limiters for upper portion of fuselage 
K interpolation index to locate wing fine grid point in global coarse 
grid 

K array to order WSI grid inner-boundary side-face J-K grid points 
K interpolation index to locate WSI inner-boundary front-face grid 
point in coarse C-grid 

K interpolation index to locate WSI inner-boundary side-face grid point 
in coarse C-grid 

K interpolation index to locate WSI grid point in global coarse grid 

Number of "coarse" grid iterations (input as AXITC) 

Number of "fine" grid iterations (input as AXITF) 

Number of "intermediate" grid iterations (input as AXITM) 

Mode of operation for wing viscous calculations (input as VISMOD) 

Number of fins in each store fin set (input as FANG) 

Fuselage input/solution flag (input as BODY) 
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NFINS 

NINB 

NINFS 

NINP 

NINS 

NINW 

NOSEB 

NOSEFS 

NOSEP 

NOSES 

NOSEW 

NPOA 

NPOS 

NPOSS 

NPYLS 

NSECFS 

NSECP 

NSECT 

NSING 

NSTOR 

NTC 

NWING 

OMEGAS 

PBL 

PCIRF 

PCIRG 

PCIRM 

PCIRS 

PCIRSS 

PCL 


Number of sets of store fins (input as FINS) 

Number of ordinates defining axi symmetric fuselage (input as BNIN) 
Number of ordinates defining each store fin section (input as FNIN) 
Number of ordinates defining each pylon section (input as PNIN) 

Number of ordinates defining store body (input as SNIN) 

Number of ordinates defining each wing section (input as ANIN) 
Blunt/sharp nose fuselage code for spline fit (input as ANOSB) 
Blunt/sharp nose store fin section code for spline fit (input as FNOSE) 
Blunt/sharp nose pylon section code for spline fit (input as PNOSE) 
Blunt/sharp nose store body code for spline fit (input as ANOSES) 
Blunt/sharp nose wing section code for spline fit (input as ANOSW) 

Number of X grid points between leading edge and trailing edge in wing 
fine grid 

Number of coarse C-grid points between store body nose and tail 
Number of fine C-grid points between store body nose and tail 
Pylon input flag (input as PYLS) 

Number of airfoils defining a set of store fins, also store fin solu- 
tion flag (input as FSEC) 

Number of defining pylon sections, also pylon solution flag (input as 
PSEC) 

Number of defining wing sections (input as ASECT) 

C-grid singularity location input flag (input as ASING) 

Store i nput/sol ution flag (input as STOR) 

Number of global coarse grid points representing fuselage cross- 
sections 

Wing input/solution flag (input as WING) 

Store quasi-steady roll rate parameter 

Wing lower surface WSI grid potentials 

Pylon circulation in wing fine grid 

Pylon circulation in global coarse grid 

Pylon circulation in WSI grid 

Pylon circulation in coarse C-grid 

Pylon circulation in fine C-grid 

Wing lower surface global coarse grid potentials 
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PFINC Fin lower surface coarse C-grid potentials 

PFINS Fin lower surface fine C-grid potentials 

PFL Wing lower surface wing fine grid potentials 

PHB WSI grid potentials 

PHC Global coarse grid potentials 

PHF Wing fine grid potentials 

PHS Coarse C-grid potentials 

PHSS Fine C-grid potentials 

PI IT 

PN0SE1 Store body nose stagnation potential in coarse C-grid 

PN0SE2 Store body nose stagnation potential in fine C-grid 

PNXC X component of pylon surface normal in coarse C-grid 

PNXF X component of pylon surface normal in wing fine grid 

PNXG X component of pylon surface normal in global coarse grid 

PNXM X component of pylon surface normal in WSI grid 

PNXS X component of pylon surface normal in fine C-grid 

PNYC R component of pylon surface normal in coarse C-grid 

PNYF Z component of pylon surface normal in wing fine grid 

PNYG Z component of pylon surface normal in global coarse grid 

PNYM Z component of pylon surface normal in WSI grid 

PNYS R component of pylon surface normal in fine C-grid 

PPYLF Pylon lower surface wing fine grid potentials 

PPYLG Pylon lower surface global coarse grid potentials 

PPYLIB Dummy wing fine grid potentials to treat pylon inboard surface 

PPYLM Pylon lower surface WSI grid potentials 

PPYLOB Dummy wing fine grid potentials to treat pylon outboard surface 

PPYLS Pylon lower surface coarse C-grid potentials 

PPYLSS Pylon lower surface fine C-grid potentials 

PRLEC R coordinate of pylon section leading edge in coarse C-grid 

PRLES R coordinate of pylon section leading edge in fine C-grid 

PRTEC R coordinate of pylon section trailing edge in coarse C-grid 

PRTES R coordinate of pylon section trailing edge in fine C-grid 

PS Pylon plot scaling parameter 

PSAVE1 Restores potentials on wing upper or lower surface plane which inter- 

polation schemes replaced with dummy values 
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PSAVE2 

PSAVE3 

PSAVE4 

PSDD 

PSFD 

PSLD 

PSMD 

PWBJ 

PWBJL 

PXLEC 

PXLEF 

PXLEG 

PXLEM 

PXLES 

PXTEC 

PXTEF 

PXTEG 

PXTEM 

PXTES 

PYCCF 

PYCCMY 

PYCCX 

PYCCZ 

PYFCF 

PYFCMY 

PYFCX 

PYFCZ 


Restores potentials above or below wing surface plane which interpola- 
tion schemes replaced with dummy values 

Restores potentials on pylon inboard or outboard surface plane which 
interpolation schemes replaced with dummy values 

Restores potentials inboard or outboard of pylon surface plane which 
interpolation schemes replaced with dummy values 

Wing spanwise drag coefficient, c^c^/e g at inboard wing/pylon junc- 
tion 

Wing spanwise skin friction coefficient, c * c f/ c ave at inboard wing/ 
pylon junction 

Wing spanwise lift coefficient, c * c ^ //c ave at inboard wing/pylon junc- 
tion 

Wing spanwise pitching moment coefficient, c * c m / c ave at inboard 
wing/pylon junction 

Dummy wing fine grid potential values inboard of wing/fuselage junction 
Dummy wing lower surface fine grid potential values inboard of 
wing/fuselage junction 

X coordinate of pylon section leading edge in coarse C-grid 

X coordinate of pylon section leading edge in wing fine grid 

X coordinate of pylon section leading edge in global coarse grid 

X coordinate of pylon section leading edge in WSI grid 

X coordinate of pylon section leading edge in fine C-grid 

X coordinate of pylon section trailing edge in coarse C-grid 

X coordinate of pylon section trailing edge in wing fine grid 

X coordinate of pylon section trailing edge in global coarse grid 

X coordinate of pylon section trailing edge in WSI grid 

X coordinate of pylon section trailing edge in fine C-grid 

Integrated pylon section skin friction coefficient in coarse C-grid 
Integrated pylon section pitching moment coefficient in coarse C-grid 
Integrated pylon section axial coefficient in coarse C-grid 
Integrated pylon section lift coefficient in coarse C-grid 
Integrated pylon section skin friction coefficient in wing fine grid 
Integrated pylon section pitching moment coefficient in wing fine grid 
Integrated pylon section axial coefficient in wing fine grid 
Integrated pylon section lift coefficient in wing fine grid 
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PYGCF 


Integrated pylon section skin friction coefficient in global coarse 
grid 

PYGCMY Integrated pylon section pitching moment coefficient in global coarse 

grid 

PYGCX Integrated pylon section axial coefficient in global coarse grid 

PYGCZ Integrated pylon section lift coefficient in global coarse grid 

PYLCMX X component of pylon moment coefficient (stability axes) 

PYLCMY Y component of pylon moment coefficient (stability axes) 

PYLCMZ Z component of pylon moment coefficient (stability axes) 

PYLCX X component of pylon force coefficient (stability axes) 

PYLCXV X component of pylon force coefficient (stability axes) including skin 

friction 

PYLCY Y component of pylon force coefficient (stability axes) 

PYLCZ Z component of pylon force coefficient (stability axes) 

PYMCF Integrated pylon section skin friction coefficient in WSI grid 

PYMCMY Integrated pylon section pitching moment coefficient in WSI grid 

PYMCX Integrated pylon section axial coefficient in WSI grid 

PYMCZ Integrated pylon section lift coefficient in WSI grid 

PYSCF Integrated pylon section skin friction coefficient in fine C-grid 

PYSCMY Integrated pylon section pitching moment coefficient in fine C-grid 

PYSCX Integrated pylon section axial coefficient in fine C-grid 

PYSCZ Integrated pylon section lift coefficient in fine C-grid 

RAV Average fuselage radius for boundary condition calculation 

RC Axisymmetric fuselage radius distribution in global coarse grid 

RE Freestream Reynolds number 

REFA Configuration reference area, used to compute force and moment coeffi- 

cients 

REFAS Store reference area, used to compute force and moment coefficients 

REFL Configuration reference length, used to compute moment coefficients 

REFLS Store reference length, used to compute moment coefficient 

REFX Configuration reference X moment center, used to compute moment coeffi- 

cients 

REFXS Store reference X moment center, used to compute moment coefficients 

RET Coarse C-grid metric r 

n 

RETET Coarse C-grid metric r 

tin 
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RIN Radial coordinates defining axisymmetric fuselage (input option) 

RINS Radial coordinates defining store body 

RLC QUICK fuselage radius distributions in global coarse grid 

RLOC Fuselage bottom centerline radius distribution 

RMAX Fuselage maximum radius 

RMAXS Store maximum radius, excluding fins 

RMAXSF Store maximum radius, including fins 

ROLLS Store roll angle 

RUPC Fuselage top centerline radius distribution 

SBDCMX X component of store body moment coefficient (stability axes) 

SBDCMY Y component of store body moment coefficient (stability axes) 

SBDCMZ Z component of store body moment coefficient (stability axes) 

SBDCX X component of store body force coefficient (stability axes) 

SBDCY Y component of store body force coefficient (stability axes) 

SBDCZ Z component of store body force coefficient (stability axes) 

SBDWET Store body wetted area 

SBVCMX X component of store body moment coefficient (stability axes) including 

skin friction and viscous crossflow 

SBVCMY Y component of store body moment coefficient (stability axes) including 

skin friction and viscous crossflow 

SBVCMZ Z component of store body moment coefficient (stability axes) including 

skin friction and viscous crossflow 

SBVCX X component of store body force coefficient (stability axes) including 

skin friction and viscous crossflow 

SBVCY Y component of store body force coefficient (stability axes) including 

skin friction and viscous crossflow 

SBVCZ Z component of store body force coefficient (stability axes) including 

skin friction and viscous crossflow 

SDD Wing spanwise drag coefficient, c * c c / c ave 

SFD Wing spanwise skin friction coefficient, c»c f /c 

T 3 V 6 

SFNWET Store fin wetted area 

SGRAD Fuselage slopes at wing/fuselage junction 

SINT Trigonometric sine of C-grid 0 coordinates 

SLD Wing spanwise lift coefficient, c»c./c 

SMD Wing spanwise pitching moment coefficient, c»c /c 

3 r r 3 m ave 
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SNOSE 

SS 

STAIL 

STRCMX 

STRCMY 

STRCMZ 

STRCX 

STRCY 

STRCZ 

STRWET 

STVCMX 

STVCMY 

STVCMZ 

STVCX 

STVCY 

STVCZ 


THETC 

THT 

TICEIO 

TICRUI 

TITLE 

TITLES 

TIWFG 

TIWSIF 


X coordinate of store body nose 
Store body plot scaling parameter 
X coordinate of store body tail 

X component of store moment coefficient (stability axes) 

Y component of store moment coefficient (stability axes) 

Z component of store moment coefficient (stability axes) 

X component of store force coefficient (stability axes) 

Y component of store force coefficient (stability axes) 

Z component of store force coefficient (stability axes) 

Store wetted area 

X component of store moment coefficient (stability axes) including skin 
friction and viscous crossflow 

Y component of store moment coefficient (stability axes) including skin 
friction and viscous crossflow 

Z component of store moment coefficient (stability axes) including skin 
friction and viscous crossflow 

X component of store force coefficient (stability axes) including skin 
friction and viscous crossflow 

Y component of store force coefficient (stability axes) including skin 
friction and viscous crossflow 

Z component of store force coefficient (stability axes) including skin 
friction and viscous crossflow 

Fuselage global coarse grid angular cuts 
C-grid 9 coordinates 

I interpolation parameter to calculate a coarse C-grid quantity based 
on WSI grid values 

I interpolation parameter to calculate a global coarse grid quantity 
based on WSI grid values 

Case title identifying print and plot output 
Store title identifying print and plot output for store 
I interpolation parameter to calculate a wing fine grid quantity based 
on global coarse grid values 

I interpolation parameter to calculate a WSI inner-boundary front-face 
grid quantity based on coarse C-grid values 
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TIWSII 


I interpolation parameter to calculate a WSI inner-boundary side-face 
grid quantity based on coarse C-grid values 
TIWSIO I interpolation parameter to calculate a WSI grid quantity based on 

global coarse grid values 

TJCE10 J interpolation parameter to calculate a coarse C-grid quantity based 

on WSI grid values 

TJCRUI J interpolation parameter to calculate a global coarse grid quantity 

based on WSI grid values 

TJWFG J interpolation parameter to calculate a wing fine grid quantity based 

on global coarse grid values 

TJWSIF J interpolation parameter to calculate a WSI inner-boundary front-face 

grid quantity based on coarse C-grid values 

TJWSII J interpolation parameter to calculate a WSI inner-boundary side-face 

grid quantity based on coarse C-grid values 
TJWSIO J interpolation parameter to calculate a WSI grid quantity based on 

global coarse grid values 

TKCE10 K interpolation parameter to calculate a coarse C-grid quantity based 

on WSI grid val ues 

TKCRUI K interpolation parameter to calculate a global coarse grid quantity 

based on WSI grid values 

TKWFG K interpolation parameter to calculate a wing fine grid quantity based 

on global coarse grid values 

TKWSIF K interpolation parameter to calculate a WSI inner-boundary front-face 

grid quantity based on coarse C-grid values 

TKWSII K interpolation parameter to calculate a WSI inner-boundary side-face 

grid quantity based on coarse C-grid values 
TKWSIO K interpolation parameter to calculate a WSI grid quantity based on 

global coarse grid values 
TSLOC Wing fine grid local sweep angles 

TWIST Wing twist (incidence) distribution 

W Relaxation factor, w 

WAREA Wing area computed from defining wing sections 

WCD Wing drag coefficient 

WCF Wing skin friction coefficient 

WCL -Wing lift coefficient 
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WCM 

WCORD 

WS 

X 

XBCS 

XBF 

XCSCUT 

XET 

XETET 

XI 
XIB 
XINB 
XINF 
XINP 
XINS 
XINW 
XIRRS 
XIRRSS 
XIRS 
XIRSS 
XIXRS 
XIXRSS 
XIXS 
XIXSS 
XIXXS 
XIXXSS 
XLE 

XLEB 

XLS 

XNC 

XOL 

XPL 

XPLF 

XPLP 


Wing pitching moment coefficient 

Global coarse and wing fine grid local wing section chord lengths 
Wing plot scaling parameter 

Global coarse grid X coordinates 

Fuselage cross-section Y coordinates (for plotting only) 

WSI grid X coordinates 

Fuselage computational cross-section Y coordinates (for plotting only) 
Coarse C-grid metric x 

n 

Coarse C-grid metric x 

tin 

Global coarse grid £ coordinates 
WSI grid £ coordinates 

X coordinates defining axi symmetric fuselage (input option) 

X ordinates defining store fin section 
X ordinates defining pylon section 
X coordinates defining store body 
X ordinates defining wing section 
Coarse C-grid metric £ 

Fine C-grid metric £ 

Coarse C-grid metric 
Fine C-grid metric £ r 
Coarse C-grid metric £ 

Fine C-grid metric £ 

Coarse C-grid metric £ x 
Fine C-grid metric £ x 
Coarse C-grid metric £ xx 
Fine C-grid metric £ 

X X 

Global coarse and wing fine grid X coordinate of local wing section 
leading edge 

WSI grid X coordinate of local wing section leading edge 
Store body length 

X component of fuselage surface normal in global coarse grid 
Non-dimensional X distance along chord or body length 
X leading edge coordinate of input wing section 
X leading edge coordinate of input store fin section 
X leading edge coordinate of input pylon section 
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XPT 

XPTF 

XPTP 

XS1 

XS2 

XSC 

XSING 

XSS 

XTE 

XTEB 

XWF 


X trailing edge coordinate of input wing section 
X trailing edge coordinate of input store fin section 
X trailing edge coordinate of input pylon section 
Upstream extent of entire wing planform 
Downstream extent of entire wing planform 
Coarse C-grid X coordinates 

X coordinate of singularity location for C-grid conformal mapping 
Fine C-grid X coordinates 

Global coarse and wing fine grid X coordinate of local wing section 
trai 1 i ng edge 

WSI grid X coordinate of local wing section trailing edge 
Wing fine grid X coordinates 


Y 

YBCS 

YBF 

YBMHB 

YCSCUT 

YFINC 

YFINS 

YINF 

YINL 

YINP 

YINU 

YNC 

YOB 

YP 

YPF 

YPYLC 

YPYLCC 

YPYLF 

YPYLG 

YPYLM 

YPYLS 

YSC 

YSFL 


Global coarse grid Y coordinates 

Fuselage cross-section Z coordinates (for plotting only) 

WSI grid Y coordinates 

Fuselage maximum half breadth coordinates (for plotting only) 

Fuselage computational cross-section Z coordinates (for plotting only) 

Y coordinate defining store fin section in coarse C-grid 

Y coordinate defining store fin section in fine C-grid 

Y ordinates defining store fin section 

Y ordinates defining lower wing section 

Y ordinates defining pylon section 

Y ordinates defining upper wing section 

Y component of fuselage surface normal in global coarse grid 
Wing span station, ti 

Y coordinate of input wing section 

Radial coordinate of input store fin section 

Y coordinate defining pylon section in coarse C-grid 

Y coordinate defining pylon section in fine C-grid 

Y coordinate defining pylon section in wing fine grid 

Y coordinate defining pylon section global coarse grid 

Y coordinate defining pylon section in WSI grid 
Spanwise Y coordinate of pylon 

Coarse C-grid R coordinates 

Y ordinates defining lower wing section at each fine grid span station 
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YSFU 

YSS 

YSTOR 

YTIP 

Z 

ZBF 

ZBODY 

ZETA 

ZETAB 

ZNC 

ZPP 

ZSTOR 

ZWF 

ZWING 


Y ordinates defining upper wing section at each fine grid span station 

Fine C-grid R coordinates 

Spanwise Y coordinate of store centerline 

Wing semi span, b/2 

Global coarse grid Z coordinates 
WSI grid Z coordinates 

Ax i symmetric fuselage Z centerline coordinate (input option) 

Global coarse grid £ coordinates 
WSI grid c coordinates 

Z component of fuselage surface normal in global coarse grid 
Z coordinate of input pylon section 
Z coordinate of store centerline 
Wing fine grid Z coordinates 
Z coordinate of wing reference plane 
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